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The Automatic Control of an Aeroplane 
in the Landing Approach onto 
an Aircraft Carrier 


J. B. HELLIWELL, Ph.D., B.Sc. 
(Department of Mathematics, College of Technology, Manchester) 


Summary: An investigation is made into the possibility of the completely 
automatic longitudinal control of an aircraft in the deck landing approach. The 
method presupposes the existence of a radio beam system transmitted from the 
carrier and a study is made of the motion of the aircraft down this beam under 
automatic control. 


The longitudinal motion alone of the aircraft and carrier is considered. The 
motion of the carrier is assumed to consist of a steady forward velocity modified by 
a non-pitching sinusoidal oscillation in the vertical plane. 


The most successful method of approach control proposed is based on the use 
of a single beam emanating from a source at a fixed distance above the carrier deck 


and fixed in direction. The preliminary pilot-controlled approach is assumed to be 
down the mean position of the beam, and the auto-control requires, over and above 
certain ideal variations of elevator angle and thrust, additional types of control of 
the following forms. On the elevator, displacement pitch and displacement track 
error control; and on the thrust, displacement track error control. 


Information is also forthcoming regarding the unsuitability of certain other 
types of control which at first sight seem promising. 


1. Introduction 


In recent years the blind landing of land-based aircraft has been given a great 
deal of attention and several systems of approach have been designed. Pritchard"? 
has given a general account of various existing types of blind approach all of which 
involve the use of radio beams and radar. The fundamental idea behind each 
different system is essentially the same, namely, to keep the aircraft on some 
predetermined flight path down onto the runway. To this end information regarding 
the position of the aircraft relative to the approach beams is periodically given to the 
pilot so that he may make certain adjustments to the controls and so centralise the 
aircraft in the beam system. 


Paper received October 1950. 
{The Aeronautical Quarterly, Volume III, February 1952] 
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The purpose of the present paper is to investigate the possibility of a similar 
technique to that just described being employed in the deck landing approach, 
In this case, however, the landing area is very limited in extent and, in general, js 
not stationary, so that a high degree of accuracy is required in the response of the 
pilot to the “error” stimuli supplied to him. Moreover, in bad weather the resources 
of the pilot will already be severely taxed and therefore automatic control is 
desirable to relieve him of the responsibility of guiding the machine down the beam, 
Furthermore, the response of an auto-pilot depends solely upon the input signals 
and thus can be assessed once and for all. From all viewpoints, therefore, it is to be 
expected that a more satisfactory approach will ensue if the “error” radio-signals 
are fed directly to the auto-pilot. The problem thus resolves into an investigation 
into the automatic control of the deck landing approach, and in this paper its 
longitudinal aspects are considered. 


The discussion separates into two distinct parts. Part I deals with the 
establishment of a beam system, and the general design of an auto-pilot suitable for 
the approach. Part II combines the general conclusions of the previous part in a 
detailed analysis of a Single Beam Approach system. 


Notation 
A initial aircraft distance from beam source 
amplitude of carrier oscillation 
¢ =[an/(V.2)] cos ¢ (see equation (1) of Appendix I) 
Cyo,Co: profile drag coefficient, flaps up and flaps down respectively 
Cy: induced drag coefficient 
Cp,Ci drag and lift coefficients of complete aircraft 
Cir _ lift coefficient of tail unit 
C mean wing chord 
e as suffix, denotes equilibrium conditions 
H aircraft height above mean deck level 
i _ as suffix, denotes ideal conditions 
iz _ coefficient of inertia about y-axis (see Ref. 4) 


k,l ideal incidence parameters (see equation (6) of Appendix I) 


ly tail arm 
k. =4C, 
N,n _ dimensional and dimensionless circular frequencies of carrier 
oscillation 


S,S’ and tailplane areas 
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GENERAL DESIGN 
BEAM LAYOUT 


AUTOMATIC 


thrust in Ib. 


time in seconds 


CONTROL 


FOR 


DECK 


LANDING 


APPROACH 


constants defining carrier displacement at auto-pilot engagement 


unit of aerodynamic time, the aerosecond 


velocity along x-axis 


increment of U from its equilibrium value 


carrier forward velocity relative to sea 


steady aircraft velocity in preliminary approach 


=(V./V.) cos » (see equation (1) of Appendix 1) 


velocity along z-axis 


aircraft weight 


increment of W from its equilibrium value 


is given by: —thrust increment in lb. = pSV .*x, 


is given by : —track error in ft. 


wing and tailplane incidences 


V.tz* 


ideal flight angle down beam (see equation (1) of Appendix I) 


flight angle down mean line of beam (see equation (2) of Appendix I) 


climbing angle of x-axis 


angle of downwash at tail 


beam angle of elevation 


increment in elevator angle 


increment in attitude 


aspect ratio 


relative aircraft density 


air density 


aerodynamic time (t=fr) 


PART I 
OF THE AUTO-PILOT AND 


2. Establishment of a Beam System 


2.1. General requirements 

First of all it is necessary to say a little about the motion of the aircraft carrier 
itself. This will generally consist of a number of oscillations about a steady state 
in which the ship steams forward on an even keel. Here all movement not confined 
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LINE OF INTERSECTION 
OF BEAMS 


————— DIRECTION OF FORWARD MOTION 


PRIMARY BEAM 
———-——- SECONDARY BEAM 
Fig. 1. 


to the longitudinal plane is neglected, and in that plane it is assumed that the carrier 
steams steadily ahead and possesses, in addition, a non-pitching sinusoidal oscillation 
in the vertical plane. Pitching as such is unimportant since only the landing 
approach is being considered and not the touch down, so that its effect on the 
attitude of the aircraft relative to the deck is of no weight and its effect on beam 
movement may be removed by suitable stabilisation. 


Now, whatever be the beam system established, the following general require- 
ments should be satisfied. First, the early stages of the approach should be a steady, 
straight-line flight with no attempt made to follow the actual oscillations of the 
carrier, and secondly, the ultimate flight path must be towards the touch-down point 
on the deck. The first stage may be under either human or automatic control; the 
second must be fully automatic. It is necessary to make sure that the ideal flight, 
in which the aircraft keeps on the beam, is practicable in the sense that excessive 
elevator movements and thrust changes or stalling do not occur. It is then further 
necessary to determine whether the auto-controlled flight remains practicable and 
stable as the machine attempts to stay on the beam in the presence of gusts and 
other disturbing phenomena. 


It might be added that throughout the complete approach the aircraft should 
be trimmed for landing conditions, namely, flaps and undercarriage down, and so on 


2.2. Twin beam system 


The fact that a distinction is made between the early and final stages of the 
approach suggests that the possibilities of a twin beam system should be investigated. 
One beam, the primary, say, should be stabilised relative to the oscillation of the 
carrier with its direction in space fixed, although still possessing the forward motion 
of the ship; this will provide a straight-line flight path for the preliminary approach. 
The other, a secondary beam, should have its source localised in the carrier deck, in 
the same vertical as that of the primary, with its direction constrained to intersect 
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AUTOMATIC CONTROI 


the primary at a constant distance, say, from the stabilised source of the latter. 
Three successive positions of the carrier and the beam system are depicted in Fig. 1. 
In the approach the pilot should maintain steady, straight-line flight down the 
primary beam and at some point of this path engage the auto-pilot. This will effect 
the change-over to the secondary beam on arrival at P and will then control the 
flight down this beam onto the deck. 


It turns out, however, that the Twin Beam Approach system must be rejected, 
as it is found that the ideal flight down the secondary beam is totally impracticable 
owing to the demand for a large negative thrust in following the seemingly 
increasing oscillations of the deck, as viewed from the aircraft. 


2.3. Single beam system 


The detailed study of the twin beam system leads to the conclusion that it 
becomes more and more satisfactory as the point of intersection of the two beams 
recedes farther away. In the limit this point is allowed to go to infinity, whereupon 
the two beams become parallel and the primary one may be dispensed with. There 
remains a single beam emanating from a source in the deck*, stabilised at a given 
inclination to the horizontal, and moving bodily up and down with the motion of the 
carrier. Such a system is fully discussed in Section 6 and an outline of the theoretical 
analysis is given in Appendix I. It will be shown that this type of blind-approach 
layout appears to be satisfactory in all essential details. 


3. Automatic Longitudinal Control about Unsteady States 
of Motion 


The automatic control of aircraft is now common practice and considerable 
literature exists on the subject. An article by Meredith and Cooke"? discusses the 
general problems of auto-control and a review of the methods employed is given 
by Weiss’. More recent papers by Neumark’ and Sudworth and Hopkin: ‘ 
considerably clarify the theoretical treatmentt. However, several important 
differences exist between the control equations developed by the foregoing and 
those required for the present work. Accordingly, the general auto-control formulae 
in any longitudinal motion will be derived with special reference to these differences. 


Analytically the longitudinal motion of an aeroplane is referred to a certain, 
steady, straight-line flight, known as Equilibrium Flight, in which an orthogonal 
system of wind axes is established. The small deviations of velocity, incidence, 
attitude and so on from their equilibrium values are known as increments. In the 


*As remarked in Section 8, the beam source does not actually need to be in the deck but, 
taken there, it simplifies the analysis considerably. 


tSince this work was concluded, Meredith®) has given a survey of miodern trends in auto- 
Pilot design, and Taylor has discussed a generalised theory of aircraft auto-control. 
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general motion one may specify that these increments should take up prescribed 
values related to corresponding settings of the controls. Thus, in the case of 
prescribed, steady, straight-line flight, the several values are taken to be zero; 
more generally they may be prescribed functions of the time. The purpose of 
auto-control is then to suppress any variations of the flighi increments from 
these ideal values. 


In what follows the dimensionless notation of Neumark’ is employed throughout 
Ideal values of the increments, together with their time rates of change, are specified 
by the suffix “i,” e.g. the attitude increment is 6. The actual values of the increments 
which occur in flight are written without the suffixes, e.g. 6. Error terms arise, 
therefore, of the type (9—4,). The control is made automatic by gearing the contro] 
movements to linear combinations of these error terms. Thus, typical formulae of 
automatic elevator and automatic thrust control become 

) d ( u )} 


dé dé, 

n=nitfo (6 —6;)+ fas AS +h( a) + fay dt { (= 
Ww 


u Uj 
Xr =Xrit Co (0-6) + ( + eu (= + 


where fo, fuer afte constant elevator control gearings and 
C6, Cab pats Cus Causats Cws Caw/at CONStant thrust control gearings. 


Automatic controls such as those devised would give the required flight with 
regard to speed and attitude but would neither be able to localise the line of flight 
in space nor to maintain it along a beam. Such a system is useless here and a 
form of control needs adding to force the aircraft to ride the beam. Hopkin’ has 
suggested a way of doing this in the case of lateral motion, introducing the concept 
of track error, and the method is now applied to the longitudinal case. 


The track error, z*, is defined as the perpendicular distance of the aircraft centre 
of gravity from some required flight path. The sign is determined by taking z 
positive in the general direction of the positive z-axis. Thus in level flight a positive 
track error means that the flight is below the required path. In the general case a 
kinematic expression can be derived relating dz*/dz to the flight increments. 


Automatic track error control is obtained by gearing the elevator and thrust 
directly to the variations of z* in the form 


n= — fas, z*_—f; 


(2) 


dz* 
X71 = Cazyat +e,2* +e: | z*dr 
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where Cts fe are additional constant thrust and elevator 
control gearings. 

Complete automatic longitudinal control is obtained by combining the track 
error controls (2) linearly with the flight controls (1). From a practical aspect such 
auto-control involves the measurement and use of recorded quantities yielding error 
terms to which the elevator and thrust are related with appropriate “ gear ratios.” 
The general relationships between these recorded quantities and the aircraft control 
settings may thus be written as 


1 
d{w w dz* 
Xr=Xni eo ( » Ome Caza €.Z e:| z*d | 


(3) 


4. Types of Auto-pilot Suitable for the Landing Approach 


The broad features of the design of the auto-pilot must be considered before an 
approach down a particular beam system is investigated in detail. The question to 
be answered is—what types of control need to be employed so that the 
“perturbation” motion, as defined in Appendix I, is stable? Owing to the labour 
involved in the calculations, only a limited number of the types of control are 
investigated; those retained may be specified by the gearings. 


It is necessary to make a simplified response calculation to obtain some idea 
of the value of any particular control gearing, and some typical values, together 
with the criteria involved in the assessment, are found to be 


dfe=22 : 0.5 deg. elevator per deg. attitude error 
df.— 16: 1 deg. elevator per 25 ft. track error 


e.=0.7 : 3 ft./sec. rate of climb per 20 ft. track error. 


4.1. Some results concerning (er, fz, fs, fas) controls 


It can be shown that in gliding flight at high incidences fz, f., faza: controls are 
always unstable. The limiting incidence for stability is given where Cy,=Coi. 
In powered flight at Jow thrusts the same result holds and the value of the critical 
incidence is still approximately that for the corresponding glide. With high thrusts, 
however, an entirely different picture may be presented. In fact the extra upward 
force due to the tilting of the thrust axis may become a dominant factor and may 
actually lead to stability (as shown later in this section). 
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It may be shown further that an e: control leads to oscillatory instability always, 
save possibly for a very narrow range of values of », bounded above and below by 
regions giving definite instability for all values of the gearing. 


4.2. Stable combinations of control 


With Low Thrust. A satisfactory solution to the problem is obtained by 
combining an fs control with an ordinary e. type. The fs control needs to be used 
to maintain a satisfactory distribution of the total damping between the various 
modes of motion, while the e, control is necessary for following the beam. It may 
be shown that, with the relevant data of Appendix II, this combination is stable for 
the following values of the gearings 


0, <1.3 dfe=20, e.<2.3 
dfe=10, e:-< 18 €:<26. 


It is easy to show that a control strength 5f,=20 is a suitable one to employ, 
and in this case the stability roots for varying e. may be tabulated as 


e. Control Subsidence Slow Mode Rapid Mode 
0 — 1.5850 + 0.4591i 1.9776 + 3.8290i 
0.5 - 2.2059 0.4756 + 0.6026 1.9840 + 3.83167 
0.75 2.3714 0.3898 + 0.81817 1.9871 + 3.83307 
1.0 — 2.4990 0.3230 + 0.96617 1.9901 + 3.83447 


Besides introducing a rapid subsidence an e. control has a marked effect on the 
frequency and damping of the slower mode. Apparently nothing is to be gained 
from the use of a strong gearing and a suitable combination should prove to be 


Se=20, e=05 . . . . . @ 


With High Thrust. The (e., fs) combination already devised is considered first. 
It may be shown that, with the data of Appendix II relevant to the present case, 
gearings 5fs=20, e.=0.75 should prove satisfactory. These values, relating to flight 
with a moderately high mean thrust, lead approximately to the same damping 
distribution as do those considered suitable for flight at low thrust. 


It has also been stated that an f, control is destabilising, except possibly in this 
case of flight with high thrust. Such control is now combined with the two fore- 
going types and, retaining the values of e.(=0.75) and df¢(=20), the stability roots 
are compared when f,=0 and 0.15. 


f. Control Subsidence Slow Mode Rapid Mode 
0 — 2.3649 — 0.4426 + 0.7958i — 1.9825 + 3.81847 
0.15 — 1.4619 — 0.8565 + 0.82887 — 2.0201 + 3.6724 
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Fig. 2. 
It is seen that the small gearing introduced has the important effect of removing 
damping from the rapid subsidence and transferring it to the slow mode. In view of 
this result these three controls should be combined in the auto-pilot with 
loy, the gearings 
7.-015,. ‘ ‘ (5) 
PART II 
A SINGLE BEAM APPROACH SYSTEM 
5. General Details 
In this part of the paper the single beam layout proposed in Section 2 and the 
auto-pilot “designed” in Section 4 are combined to give a single beam approach 
. system. The scheme envisaged consists of a single beam, fixed in direction, making a 


constant angle » with the carrier deck. Its source is localised in the deck, thereby 
(4) possessing the complete carrier motion. Fig. 2 shows three successive positions of 
the carrier and beam during the motion. 


At some large distance from the carrier the aircraft should be trimmed for 


wr the landing approach and, under the control of the pilot, should be guided towards 
i the beam and maintained in flight down the mean line of it. At some stage of this 


flight the auto-pilot should be engaged and this will then direct the aircraft into the 
oscillating beam and down the beam, on to the deck. The theoretical analysis of 
this the complete flight with this preliminary approach is outlined in Appendix I. 


Sf 6. Ideal Flight Down The Beam 


The formal expressions for the ideal increments are given in Appendix I as 
the set of equations (10)-(14). The variations of flight from that down the mean 
line of the beam are of a sinusoidal nature of the same frequency as the carrier 
oscillation. Consider the thrust. It is clearly important to maintain a large thrust 
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NUMERALS 1,2 AND 3 INDICATE CORRESPONDING POSITIONS OF AIRCRAFT AND TOUCH=DOWN 
POINT AT SUCCESSIVE INTERVALS 


SECTION E-F OF THE IDEAL FLIGHT PATH 


AIRCRAFT AT 


F E 


TOUCH-DOWN AT 
87-5 SECS. 


3700 FT. > 2080 - 
TOUCH: DOWN POINT AT #0 


PATH OF AIRCRAFT PATH OF TOUCH-DOWN POINT 
Fig. 3. 


in the preliminary approach down the mean line of the beam since the amplitude 
of the ideal thrust oscillation may be of considerable magnitude, and a large negative 
value is incapable of being produced by any conventional aircraft. This implies 
that the beam angle must be small and the mean angle of flight, 8,,, must be 
considerably less than the gliding angle. Such a value is taken in the numerical 
data of Appendix II. 


Further, it may be shown that it is necessary to incorporate a prescribed 
variation of the incidence in the ideal flight. With a sinusoidally oscillating incidence 
with parameters k and / as laid down in Appendix I, the amplitudes of the ideal 
flight increments are themselves functions of k and /. These parameters are chosen 
so as to specify such incidence changes as cause the smallest corresponding 
oscillations of speed, attitude and control movements consistent with not stalling 
the aircraft as it flies down the beam. 


Values of k and / which satisfy these requirements are 
k=-0.06 /=0.01 . 6 


and the ideal flight down the beam takes place with the following oscillation of 
the type dcos(nt+¢) and period 8 sec. 


Amplitude (d) Mean Setting 


Incidence: deg. 3.5 12 
Speed: ft./sec. 3.9 123.6 
Attitude: deg. 5.4 —4.9 
Thrust: Ib. 370 1090 


Elevator angle: deg. 2.2 — 


y M 
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Moreover, defining the normal acceleration as the reversed acceleration along the 
z-axis, the maximum excess value caused by the manoeuvres is about 0.2 g. 


Thus an ideal flight down the beam is obtained, capable of practical achieve- 
ment by a conventional aircraft. The flight path relative to the touch-down point on 
the deck is drawn in Fig. 3 and the variations of the flight increments and control 
settings are plotted in Figs. 4 to 8. 


7. Actual Flight Down the Beam with (e:, fs, f-) Auto-control 


Having determined a satisfactory ideal motion for the aeroplane there remains 
the final design of an automatic control system which will guide the aircraft into 
the beam and cause its motion actually to approach the ideal. Alternative stable 
combinations of control have been examined in Section 4 and their application will 
now be discussed. Here the ideal flight demands a high mean thrust and thus a 
combination incorporating f, control may be used and so the beam-riding element 
may be strengthened. Therefore an auto-pilot combining (e., fs, f.) controls with the 
stable set of gear ratios determined previously (equation (5)) is considered. 


The actual flight is now calculated as the sum of the ideal motion and the 
transient perturbations. The ideal motion has been obtained already and discussed 
in Section 6. The transient perturbations are calculated as the complete solutions 
of equations (14) of the Appendix. 

Three sets of non-ideal initial conditions at auto-pilot engagement are taken 
in the assessment of the perturbations, and in all cases the unfavourable initial 
displacement of the carrier deck from its mean level is assumed to be a maximum. 


83 
\ 
/ 
7s 


= 3 
w 
SY 


FLIGHT ACTUAL FLIGHT 4CASE (ii) ——-———~ 
casei) 


Ide 
ive 
lies 

be 
cal 
ed 
ice 
eal 
en 
ing 
KNOTS } 
(6) 
of 
Fig. 4. 
251 


J. B. HELLIWELL 


— 


(I) ------- 
IDEAL FLIGHT: ACTUAL FLIGHT < CASE (ii) —— —— 
CASE (iii) —-——.- 
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Fig. 9. 


Of the cases chosen, the first is the best that could be expected to occur in practice; 
the remainder introduce serious flight errors at auto-pilot engagement. They are 


(i) Aircraft in equilibrium flight down the mean line of the beam. 


(ii) Aircraft in flight down the mean line of the beam with a speed error from 
equilibrium of 10 knots (actual speed = 83 knots). 


(iii) Aircraft in equilibrium flight parallel to that down the mean beam with 
a 50 ft. displacement error (track error=50 ft. + initial beam 
displacement). 


The results of the calculations of the perturbations, and thus of the actual 
motion, in each of these cases is plotted in Figs. 4 to 8. The track error is plotted 
separately in Fig. 9, while Fig. 10 shows the actual flight path of the aircraft down 
the beam in the ideal case, together with the two worst cases of initial conditions, 
plotted until the perturbations have disappeared. 


A study of the curves leads to the following conclusions. In all cases a 
satisfactory response ensues, the stall is avoided and ample thrust is available. The 
total time of flight from auto-pilot engagement at about 1,000 ft. to touch-down on 
the deck is approximately 87.5 sec., and for the last 55 seconds of this the approach 
takes place ideally. In practice, conditions at auto-pilot engagement will generally 
be better than those discussed here and in such cases the approach will certainly 
be possible. 
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On an actual aircraft the presence of friction and inertia in the control circuits 
will remove the “kicks” that appear in a theoretical treatment of displacement 
type controls. It may be shown that a viscous friction type of control lag does not 
affect the foregoing conclusions adversely; in fact it merely removes these kicks 
from the theoretical response. 


Finally, as may be expected, it can be shown that the effect of a headwind is to 
bring about an all-round improvement in the approach. 


8. Some Practical Aspects 

The paper as a whole presents a description of a possible approach technique, 
with but a single example worked out, rather than a final design for an auto- 
controlled approach system. For example, it may be considered that the total time 
of approach-flight, namely 87.5 sec., is too long for practical purposes. This defect 
may be remedied by slightly modifying the approach servos so that, for example, 
the auto-pilot engagement takes place at 500 ft. instead of at 1,000 ft. as in the 
case worked out. This would result in a total approach time of about 44 sec. 
with the last 11 sec. ideal. 


Again, the auto-pilot gearings as finally given numerically in the paper are 
certainly not the best suited to the purpose. Taylor'’’ has shown that the “best” 
response to any input signal is produced when the gearings are designed for critical 
damping—best in the sense of least peak load or overshoot for a given transient 
decay time—and clearly such a modification would help in the present problem. A 
possible result would be to make it possible to lower the auto-pilot engagement 
height still further, thus obtaining a smaller total approach-flight time retaining the 
same final 11 seconds of idealised approach. 


The analysis is based upon the fact that a prescribed variation of incidence in 
ideal flight down the beam demands a certain corresponding thrust and elevator 
setting. It may be shown that, for any particular aircraft, these resulting control 
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movements may be written purely as functions of the carrier oscillation so long 
as the violence of this is less than that specified in Appendix II. The basic servo. 
mechanisms are therefore relatively simple, being merely sinusoidal functions of 
two variables. 


Next, although the beam source has been taken in the deck it is unnecessary 
to persist in such an awkward location. A suitable position would be at the 
masthead and this has the added advantage of cutting down “ground return.” The 
track error control would then be so designed as to centre on a line parallel to 
the guiding beam. 


In conclusion it is remarked that although a simple carrier oscillation has been 
assumed, the approach system is applicable to more general carrier motions provided 
that these can always be resolved by Fourier Analysis into sinusoidal non-pitching 
oscillations in the vertical plane. 
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APPENDIX I 


Theoretical Analysis of a Single Beam Approach System with a 
Preliminary Approach Down the Mean Line of the Beam 


The beam approach system and the approach flight to be analysed are 
described in Section 5 and the assumed carrier motion is as stated in Section 2 of the 
main paper. A first approximation is taken throughout the analysis. 

Two distinct problems are involved : — 

To determine the nature of the motion down the beam perfectly executed. 
To find what is the motion possessed by the aircraft as the auto-pilot controls it 
in an attempt to follow the beam. 

The following definitions are assumed : — 


Ideal Motion consists of the steady flight down the mean line of the beam and the 
subsequent perfect motion with a prescribed variation of incidence as the aircraft 
flies down the oscillating beam. 


Actual Motion consists of the auto-pilot engagement and the subsequent flight as 
the aircraft attempts to follow the oscillating beam. 
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1. Kinematic Preliminaries 


On the assumption that the aircraft c.g. remains on the beam at all times in 
ideal flight, the angle of descent down the beam may be determined from purely 
geometrical considerations. Fig. A.1 shows the positions of the beam source 
(D, D’), beam (DS, D’S’) and aircraft (S, S’) at two consecutive instants of time, with 
the vertical carrier displacement at any instant expressed as asin N (t+1,). 


aSIN N(t+§t+t,) a SIN N(t+t,) 


Fig. A.1. 


Thus, by simple geometry, if (V. sing) and (aN) are small compared with V., and 
is a small angle, the downward angle of flight along the oscillating beam becomes 


B= o—V ccosn(r+7,), . (1) 

the downward angle of flight along the mean line of the beam is 
Bu=o V tan®@, é : (2) 

and the height above mean carrier deck level may be deduced as 
 @) 


2. Axes of Reference 


Equilibrium wind axes and velocities are chosen at any convenient instant 
during the steady straight-line preliminary approach. Fig. A.2 shows the system 
of axes obtained. 


The notation is standard (Ref. 4) and we have 


y= — Ba o+V tang, 
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TAN 


Fig. A.2. to | 


3. Conditions of Ideal Flight 


The ideal flight can be specified in terms of the velocity and attitude increments, tota 
(u,w,9), from the equilibrium flight just defined. Then the elevator setting and 
thrust may be determined so that the aircraft may perform this required motion. 

The longitudinal equations of motion of an aeroplane are used here in assessing these 
ideal increments and, following Neumark“, are written 


d ( u ) u ) w ) 6b 
U. Xe +k, (4) (xr)=0 
i+) dr U. <w U. = dz +k, tan Ye (9) = 0 33 
Net 
(5) 
Of the five increments occurring in these three equations two may apparently be i 
chosen arbitrarily. One of these, however, is fixed by the previous specification that 
the aircraft c.g. remains on the beam thereby relating #, and w,, and for ease of 
analysis the other is taken to be a prescribed variation of incidence. 34 
3.1. Variation of the incidence subs 


The beam oscillates with a frequency n. It is assumed, therefore, that in flying 
ideally down the beam the incidence is prescribed to oscillate about its steady, 
preliminary approach value, with the same frequency, so that the small incidence 
variations can be written 
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w 
Then, since U, incidence increment, 
e 


the expression for w,/U, is determined immediately from equation (6) as 


3.2. Variation of the attitude 


The motion of the aircraft c.g. is such that it remains on the oscillating beam 
at all times. The constant, preliminary approach c.g. velocity vector has been taken 
to define the location of the x-axis in the body and its equilibrium direction. 
Assuming for the moment no variation of the incidence, the x-axis remains along 
the c.g. velocity vector, which is now variable, and a variation of attitude ensues, 
4, say. A further contribution, 6, say, arises from the incidence changes so that the 
total attitude change is + 


Now 4 is measured positively for an increase in climbing angle and thus 
6,= —(8-8,)=ccosn(r+7,). 
4, is simply equal to the incidence change so that the complete expression for 
4, becomes 


6,=ksinn(z7+7,)+ (1+ c)cosn(++7,). (8) 


3.3. Variation of the velocity 


Having now found expressions for w, and 4,, the second equation of motion 
(equations (5)) is used to determine u;. Introducing the portmanteau notation of 
Neumark"? this is 


uy 


{ [K(N, —Q,)+en] sin n(>+7,)+[1(N, —Q,)+ck, tan y.] cosn(7+7,) } 


(9) 


34. Variations of thrust and elevator settings 


The ideal thrust and elevator increments x;; and 1, are obtained directly by the 


| substitution of u;, wi, 4 in the first and third of equations (5). They are 


Xn= Like, 
. (10) 


[k (n? + In (v +x) + env] sin n (= + 7,) 


-[kn(¥ + \)—1(n? — cos n(= +7,) } ‘ 
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) sin 
COS n ) 
(T+1, do' 
du* 
du" 
dr 
e 
Fig. A.3. 
4. Relation Between Track Error and Flight Increments 
At any instant the position and dimensionless velocity of the aircraft c.g. relative 
to the beam are as shown in Fig. A.3, where G is the actual c.g. position and G Th 
is its projection upon the beam, so that G,G is the track error z* (now in mo 
dimensionless units) the 
Then dz*/dz is the dimensionless velocity of G relative to G, in direction G,G, con 
which may be shown to be dire 
and 
1z* ’ 
5. Equations of Controlled Motion Down the Beam Air 
The controlled equations of motion are obtained by combining the equations 
of auto-control (equation (3) of the main paper and equation (12) of this Appendix) 
with the three equations of motion (5). These six equations may be simplified by the 
use of a certain transformation relating the actual and ideal motions and defining 
their difference as “perturbations.” Thus, denoting these by * in the following 
scheme, we write 
Ler 
+u*U, 6=6, +0* | 
w=w,+wl, (13) 
+0 
It is clear that the ideal motion satisfies the equations of motion (5) and thus Car 
a simple substitution of equations (13) leads to the same set (5) with u* replacing y - 
; u/U., and so on. Similarly the control equations (equations (3) of the main paper) 
simplify into relations between the perturbations, and likewise equation (12) becomes 
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explicitly independent of the ideal motion. A slight rearrangement of these 
simplified equations leads to the following set giving the controlled motion 
down the beam. 


du* dw* d6* dz* 
+ e,u™ Caw dt dr_ e,w* + dt Caz dt dr T + ez z*dr =0 
7 7 aT 
du* 
+k.0* -xz* =0 
dr 
dw* 
-Zu*+ —Z,w* - - +k, tan y.. =@ 
dr dr 
dw* d*6* 
dr dz? dr 
du* , dw do* dz* 
— — faw sar — + n* + | z*dr 
dz dr dr dr 
— yw? + dz* =0 
dr 
(14) 


These “ perturbation” equations now contain no explicit reference to the ideal 
motion and thus the design of a satisfactory auto-pilot may be entirely removed from 
the remainder of the problem and treated as a separate entity. Any resulting 
conclusions would be of general validity; for example, the results could be applied 
directly to the case of an airfield blind approach for there the beam is stationary 
and no ideal motion exists. 


APPENDIX Il 


Numerical Data 
Aircraft Design Details 


Wing loading=40 Ib./ft.*, 1/29 ig= 1/12. 
S=250 ft.’, 0a 
aC, A=6, =0.4. 
~——=45, 
da S’/S=1/6, 
c=6 ft., OCir 5 
Length of tail arm /;=18 ft., 
C,=C, (flaps up)+ AC, (flaps) =10 +12 =22. 
Cyo=Cor + Coi =0.15+0.28 =0.43. 
W a W 
Carrier and Beam System 
V.=25 knots, a=10 ft., N =8 sec., @=15 deg. (low thrust). 


@==7.5 deg. (high thrust). 
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Equilibrium Flight Conditions (Preliminary Approach) 
V.=73 knots., a=12 deg., Y= —9.8 deg., T.=233 Ib. (low thrust), 
ye= —4.9 deg., T.= 1090 Ib. (high thrust), 
Aerodynamic Derivatives 
o=2, v=3, x=1.2, 5=43:5. 
X,= — 0.46 (low thrust), x,=—0.55 high thrust). 
X0=0.52, w=-2.2, z= -—2.465, k.=1.1. 


Other Details 

Height above mean sea level at auto-pilot engagement =1,000 ft. 

Mean height of deck above mean sea level = 60 ft. 

Ideal time of flight down the beam =87.5 sec. 
N,, P;, Q,, R,, T, portmanteau symbols (see Ref. 4) 
thrust control gearings 

: See equations (1) and (2) 

elevator control gearings 
wind axes of reference, fixed in body 
Kus Xu» Zan Ze aerodynamic force derivatives (see Ref. 4) 
aerodynamic modified moment derivatives (see Ref. 4) 


u*, w*, 6*, x;*, »* perturbation variables (see equation (13) of Appendix) 
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Supersonic Flow over Thin Symmetrical 
Wings with Given Surface Pressure 
Distribution 


F. A. GOLDSWORTHY 
(Department of Mathematics, The University, Manchester) 


Summary: The inverse problem of determining the supersonic flow past a thin 
symmetrical wing at zero incidence with given surface pressure distribution is solved 
for wings of arbitrary plan form. Expressions are obtained for the vertical derivative 
of the potential on the wing surface from which, using the linearised boundary 
condition of zero normal velocity at the surface, the profile of the wing can be 
designed. The integral equation method adopted by J. C. Evvard and extended by 
G.N. Ward is used. The analysis cannot be applied to pointed wings, whose leading 
edges are subsonic. The results in Part I are completely general and are applied to 
specific problems in Part II. 


|. Introduction 


The inverse problem of designing the profile of a thin wing set at zero incidence 
to an initially uniform stream and with given pressure distribution on the surface 
is the last of a set of four problems to be solved in the linearised theory of supersonic 
flow past thin wings. The two direct problems of determining the flow past a thin 
wing of given profile have been treated by several authors, in particular by Puckett"? 


_ and Evvard’’ for symmetrical flow, and by Ward‘ in the anti-symmetrical case 


(i.e, flow at incidence over flat plates). The determination of the flow in the former 


| of these two problems is easy in principle since the boundary conditions over the 


whole plane, in which the wing lies, are not mixed. The same is also true for the 
inverse problem in anti-symmetrical flow. The inverse problem for the symmetrical 


The gas is assumed perfect and the flow taken to be inviscid, irrotational, and 


} isentropic; then, neglecting the squares and higher powers of the disturbance 


Velocities due to the presence of the wing, linearised theory may be applied. The 


) tesulting linearised equation of motion is solved by using a result due to Hadamard’, 


and the wing profile is designed by using the linearised boundary condition at the 
wing surface. 


Paper received June 1950. 


|The Aeronautical Quarterly, Volume III, Febreary 1952] 
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The method of solution is similar to that used by Ward, integral equations being 
established to determine the velocity normal to the plane of the wing in terms of 
the given pressure distribution on the wing surface. Off the wing the normal velocity 
is known and the pressure distribution may be determined by solving further 
integral equations. By using the transformation to characteristic co-ordinates due 
to Evvard these expressions are written as integral equations of the Abel type 
and therefore admit ease in solving. 


Part II of this paper is concerned with the application of the general theory to 
two specific problems of designing the profiles of wings of rectangular and arrowhead 
plan forms, on whose surface linear pressure distributions are taken. Singularities 
enter into the solution of both problems and, in the case of the rectangular wing, 
can be explained by taking the tip edges to be rounded. The result for the arrowhead 
wing clearly demonstrates that discontinuous pressure gradients give rise to 
singular solutions. 


Notation 
x,y,z Cartesian co-ordinates with the x-axis parallel to the direction 
of the main stream 
U _ velocity of the main stream 
p, density of the main stream 
M Mach number of the main stream 
B_ defined by B?=M?-1 
u_— disturbance velocity due to the presence of the wing 
@ velocity potential defined by equation (2) 
y co-ordinates of the points in the area of integration given 
by inequality (4) 
z=Z(x,y) equation of the wing surface 
p (Part I) disturbed pressure over the wing surface 
p (Part II) Heaviside operator 
2,8 characteristic co-ordinates given by z=x— By, B=x-+ By 
a’, 8’ characteristic co-ordinates given by a’=x’— By’, B’=x'+ By 
L(2,8) defined as L(2,8)= —(2 B/x)¢ (x, y,+0) 
S projected area of the wing onto the plane z=0 
C boundary of area S 
R__ remainder of plane z=0 
O origin of co-ordinates 
a,,8, point of intersection P (see Fig. 1) 
A, (z),B,(z) functions of a defined by equation (15) 
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A, (8), B,(8) functions of defined by equation (16) 
f(z,8) the value of L(2, 8) on § 


2:,8.,h,,h, the values of L (2,8) on special regions of R (see equations 
(18) and (19)) 


F, (z),F,(z) functions introduced in equation (22) 


parameter introduced in solving equation (23) 


ret 


£,£’,n’ variables defined by equation (41) 
x=a_ equation of leading edge of the rectangular wing 
wing chord 

~=breadth of wing 
: constants in equation (37) 
constants in equation (38) 
constants in equation (49) 

7 angle of sweepback (see Fig. 5) 
fe.f. values of f(2,8) on the right side and the left side of the 

arrowhead wing respectively 


=, regions defined in Fig. 5 


PART I 
GENERAL THEORY 


2. Fundamental Equations and Boundary Conditions 

A wing, whose surface makes everywhere a small angle with the x- and y-axes 
so that it lies approximately in the plane z=0 in suitable Cartesian co-ordinates, 
is placed in a uniform supersonic stream parallel to the x-axis with velocity U, 
density p, and Mach number M. Neglecting viscosity and heat conduction the flow 
is practically isentropic (since the change of entropy at weak shocks is proportional 
to the cube of their strength) so the linearised equation of isentropic flow 

- + «4s = 
dx? dy? dz? 

holds as a first approximation, where B* = M* — 1 and the disturbance velocity to the 
uniform stream is given by 


Equation (1) is a second order hyperbolic equation and the solution can be 
given in terms of data given on the plane z=0 in a case when u=0 for x less than 
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some constant. Assuming that the potential on z= +0, 9(x, y,0), is continuoys 
and bounded, the solution of equation (1) for z>0, say, is 


where the integration is taken over the area of the plane z=0 for which 
Expression (3) is obtained by differentiating with respect to z the equation for 


{ » dz in terms of the values of ¢ on z= +0. This latter equation is obtained by 
0 z 


substitution of dz for in Hadamard’s well-known result'*’ expressing in 


0 
terms of the values of 0¢/0z on z= +0. 


If the wing is symmetric, the disturbance potential and its x- and y-derivatives 
are even functions of z, and 09/0z is an odd function of z. Let the equation of the 
wing surface be z=Z (x,y), then the boundary condition of zero normal velocity at 
the upper surface, to a linear approximation, is 


On the rest of the plane z=0, 09/0z, as a continuous odd function of z, must be zero 


When the disturbed pressure distribution, p, is specified over the wing surface, 
as it will be in this case, integration with respect to x of the linearised Bernoulli 
equation 


determines the value of 9 (x, y,+0) over the wing surface. 


Summarising, the disturbance potential 9 (x, y, z) satisfies equation (1) and in 
the region z>0, say, has the following mixed boundary conditions: 


(i) ¢ is continuous and bounded everywhere in z>+0; 
(ii) ¢>=0 ahead of some plane x=constant; 
(iii) # (x, y, +0) takes given values on the wing; 

(iv) d¢/8z=0 on the plane z=0, except over the wing. 


Then, if 0¢/0z is found over the wing, using condition (5) the wing can be designed 
to produce the required pressure distribution. 
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éo/0z can best be calculated by differentiating equation (3) with respect to z: — 


7 


Since the quantity contained in the square brackets is the integral of ¢ with respect 
to z and therefore satisfies equation (1), equation (8) may be written 


= __ +0) dx’ dy” ] 


Now taking the special case of equation (9) when z= +0, the value of 0¢/@z on 
z= +0 is given by 


(<2), +0 Ox? dy’ J J (10) 


the integration being taken over the area of the plane z=0 for which 


The form of expression (10) is appropriate for the application of the 
transformation into characteristic co-ordinates (z, 8) given by 


Putting (2 B/=)o(x’,y’,+0)= —L(2’,8’) and making the transformation (12), 
equation (10) becomes 


L (2’, da’ dp’ 
0208 LJ { (2-2')(8- BY} J (13) 
the limits of integration being 


In the manipulation of equation (13) it is convenient to divide the wing surface 
as shown in Fig. 1. 


Let the projection of the wing on to the plane z=0 be S, with boundary C, and 
denote the remainder of the plane z=0 by R. In Fig. 1, L, P,, L,P,, T,U, T,U 
are characteristic lines, « or 8=constant, for equation (1) in the (x,y) plane and 
tangents to the curve C at points L,, L,, T,, T, respectively. L, P, L, P,T,U,,T,U, 
are also characteristic lines. The origin of co-ordinates is taken at O, the point of 
intersection of L,P, and L,P,, the co-ordinates of L,, L, being (2, 0), (0, 8,) 
respectively. The curvature of C is assumed to have the same sign everywhere 
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Fig. 1. 


(i.e. C is convex), so that it may be defined in terms of single-valued functions : — 
On T,L,L,, B=A,(@): on L,7,7,. B=A,(e),. . & 
or alternatively : 
On L, L.T., (8); on T,7,L,, z=B,(8). . (16) 


Integral equations are now set up to determine 09¢/0z on S, this being done 
in intermediate stages. 


3. The General Wing 


A general case, as indicated in Fig. 2, is considered immediately. 

The problem is fully solved without real difficulty when the wing has aspect ratio 
sufficiently large for the subsonic edges L, 7, and L,T, to be independent in the 
sense that the characteristic lines L,P, L,P intersect C behind them, i.e. in the 
supersonic trailing edge T,7T,. A special case of this is when the wing-tip flows 
are independent, i.e. when the characteristic lines L,P, L,P intersect at a point 
P lying outside S. If the aspect ratio is too small the flow is theoretically 
determinable, but difficulty in solving the integral equations is encountered. 
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The range of influence of the disturbance due to the wing on the plane z= 
is the region downstream from P,L,L,P,. Elsewhere on the plane z=0 the 
potential vanishes identically. Hence integration in equation (13) need only be 
taken over that part of the region given by inequalities (14) which lies downstream 


from P, L, L, 
For the wing plan form shown in Fig. 2, let 
L(2z,8)=f (2,8) on S; 
=g,(z,8) in region P, L, M,Q,; 
=g,(z,8) in region 
=h, (2,8) in region 
=h, (2,8) in region 0,M,N.R., 


where functions g, (2, 8), h, (2, 8), and so on are, as yet, unknown. 
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(i) For a point (z,8) on S in region a (see Fig. 2) the z-derivative of the potential 
is obtained from equation (13) by inserting limits derived from (15) and (16) and is 


B,(8) Aa’) 


(ii) For a point in region P,L,M,Q. (i.e. Wa<2,; B>A,(2z)) we have from 
equation (13) and condition (7) (iv) 


Aa’) 


Hence 


(B—-B) 


A A fa’) 


for some functions F, (2), F, (8). 


Equation (22) is an integral equation of the Abel type for the quantity contained 
in the curly brackets. By condition (7) (i) f and g, are bounded; hence the expression 
in the curly brackets is bounded and F, (0)+F,(8)=0. The solution under these 
conditions is 


A fa) B « 

| | | F,’ (2’) da’ 

0 


A,(2) 
Letting 8—>A, (z), the right hand side is determined as 
A fa) 


Equation (23), together with result (24), forms an Abel integral equation for 
g. (2, 8) with solution 


a) 


(25) 


| 
1 
( 
| 
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From equation (25) the final expression for g, (2,8) is determined as 


f (2, 8’) dB’ 


(26) 


82 (2, B)= 


It can easily be verified that the potential is continuous at 8 = A, (2) in the sense that 


lim g, (2, 8)=f (2, A, (2)). 
BA fa) 
Thus the value of the potential in region P, L,M,Q, is given by expression 
(26) and is continuous at the edges of the wing. Similarly the potential at a point 
in the region P, L, M, Q, (see Fig. 2) can be determined. 


(iii) Using expressions (13), 0¢/0z at a point (2, 8) in region a, of S (i.e. 0<2< 2; 
3,<B<A,(2)) is given by 


a B 
dz/--40 O20BL) V(B-B) 
B,(8) A f(a’) 
f (2", By dB" | 82 28 
A(x’) A,(2’) 


Applying equations (23) and (24), the latter term of this equation can be 
replaced by a term involving the known function, f(z, 8), only; thus equation (28) 
becomes 


a Bg) A,fa’) 
0z/:-40 0208 V(a-2) V(B-B) V(a-2’) ) { Az 
Bs) Aa’) A fa’) 
(29) 
or at points where 0B,(8)/08 is continuous, 
(8) A,(a’) 


since A,{B,(8)}=8. 


A similar expression can be obtained for a point in region a, (Fig. 2), namely 


B a 


A (a) 
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When the flows over the wing-tips are independent and the pressure distribution 
on the wing surface is given, then the z-derivative of the potential can be calculated 
everywhere on S using expressions (20), (30), and (31), and the profile of the wing 
determined. If this is not the case then further calculations are necessary, as 
follows : — 


(iv) For a point (a, 8) on S in region b expression (13) becomes : — 


A (2) 


a 
) 


A (a) B,{A,(2)} A,(2’) 


Bg) Aa’) 


A, (2) Aja’) 


A (2) Bs’) 


0 B(8) 


Of the integrals in equation (32) the first is taken over the lower shaded rectangle 
in Fig. 3, the second over the upper shaded rectangle; the others cover the remainder 
of the dependence domain, and each also contains the upper shaded rectangle. The 
third integral occurs in equation (28) and was simplified to become the second term 
in equation (29); similarly the fourth integral can be simplified. Using these results, 
expression (32) becomes :— 


8 B(g) 
Oz J (2—2’) (B-B’) J (2—2’) J 
B,(8) A,() B,{A,(2)} A,(e’) 
0? f (B, (9), 
A,{B\p)} 
{ B—A, (2) } ) ( A, (a) (33) 
B,{A,@)} 


Thus (09/0z)._,9 is given by expression (13) taken over the lower shaded 
area, and (with the sign changed) over the upper shaded area, plus two line integrals 
along the characteristic lines which intersect at the point K (which has co-ordinates 
B,(8), A,(z)) in region @ and bound the areas. 


The first two terms of expression (33) are analogous to the result obtained 
by Ward? for the potential at a point on a wing of zero thickness and given shape. 
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It is evident from Fig. 2 that, if the flows round the subsonic edges are 
independent, regions a, @,, @,, b cover the whole wing; thus expression (33) and 
the preceeding ones determine the value of 0¢/0z everywhere on the wing surface. 
No case where this is not so will be worked out. An indication is given in (v), 
which follows, how the calculation could be progressively continued in such cases 
by methods analogous to those that have gone before. 


(v) For a point (2, 8) in region Q, M, N, R, (Fig. 2) we have from equation (13) and 
condition (7) (iv) 


z A Ja’) B 


d20BL V(a—a’) V(B-B) V(B-B) J * 
0 Aa’) A 
A,(a’) B 
+ _ dat 8: (2, dB" f (2, dp" | |: 
a, 0 A,(2’) 
(34) 


The first term of this equation is zero, for the expression in the curly brackets, 
by equations (23) and (24), is independent of 8. An Abel integral equation remains 
from which, using the method of (ii) of this Section, the potential h, (a, 8) in the 
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region under consideration can be found. The resultant expression is 


A,(e) A (a) 
_ B-A,(2)} f (2, B’) dB’ g, (2’, B’) dB’ 
h; (2, B)= { | A, (@)—B’}* A, (2)— } 


A ,(a) 0 


Similarly the potential in region Q, M, N, R, can be determined and is 
B, (8) } f(a’, B)d (2’, B’)d 
_ vi 2-B,(§)} 82 (2’, B’) da’ 


B, 
(36) 


B,(g) 


Since g, (2, 8), g, (a, 8) are known functions from equation (26) and the corresponding 
equation for the left wing tip, /, (z, 8), h, (z, 8) are determined from equations (35) 
and (36). 


By application of equation (13) the value of 09/0z for a point in any one of 
the regions 5,,b, and c (Fig. 2) can be expressed in terms of known functions. 


If N, or N, lies outside T,7,, the regions of the wing lettered in Fig. 2 do 
not cover it completely and the procedure indicated in (ii), (iii) and (iv) of this 
Section has to be repeated yet again, systematically finding the value of ¢ off the 
wing, then 0¢/0z on the wing surface. It is seen that for wings of very small aspect 
ratio the method will be extremely laborious. 


In the limiting case when the wing is pointed and the leading edges lie inside 
the Mach cone from the apex (i.e. are “ subsonic”) the points L, and L, coincide and 
the method cannot be applied at all. 


PART Il 
SPECIFIC PROBLEMS 


4. The Rectangular Wing 


In this section a wing of rectangular plan form, with leading edge normal to 
the direction of the undisturbed flow, is considered; the isobars are taken parallel 
to the leading edge and along any vertical section of the wing the pressure varies 
linearly. 


Let the wing have chord of length d and breadth 2/ such that the ratio (d//)< B 
(i.e. independent wing-tip flows). Taking the isobars to be lines x=constant and 
using equation (6) together with the given linear pressure law, the potential at the 
wing surface is found to be of the form 
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where C,, D, are known constants and f(y) is a function determined from the 
condition that the potential vanishes at the leading edge x=a. 


The final expression for the potential on the wing surface becomes 
or in characteristic co-ordinates, 
@=(C,/4)(z+ 8 . (39) 
where C, and D, are known constants. 
Using expression (20), d¢/0z at a point (2, 8) in region L, M, M, L, (Fig. 4), is 


B 
B dz’ [((C,/2) (2+ 8’ —2a)+D, (2+ 8’ —2a)] dB"), 


Vv (8-8) 


(40) 
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By substituting new variables given by 
=a’ + B’—2a, &’=2'+B-2a, €=2+B-2a, 


expression (40) reduces to: — 


(2 
0z 


which is an integral of the Faltung type, so that Heaviside operational methods 
may be applied : — 


= —~B{2C,(x-a)+D,}. 


42) 
VE 


oe 
0g 
0 


Thus the equation of the wing surface, using condition (5), is given by 
z=Z (x, y)= —B { C,x*-2C,ax+D,x }'+2(y), ; (46) 


where g(y) is a function to be determined from the condition that the wing closes 
up at the leading edge. The wing must also close up at the trailing edge thus 
imposing a condition between the constants C, and D,. Using these results 
expression (46) becomes 


z=BC, (x-a)(a+d-x), . (47) 


which checks with ordinary two-dimensional theory. 


For a point in L,M,T, the value of 0¢/0z obtained from expression (30) is 


=| { 2C,(x-a)+D, } /{ 


(48) 


As the wing edge L, T, is approached, it is seen that 0¢/0z behaves singularly like 
the inverse square root of the distance from the edge. This implies that the wing 
surface has a vertical tangent there (i.e. the tip edges are rounded). The linearised 
theory, as given in Part I of this paper, cannot be expected to hold at the rounded 
edges; but the implication that the tip edges are rounded is perhaps more than 
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could have been expected from a linearised theory, for it is in good agreement 
with the design requirements of subsonic edges. Thus, from equation (5) the wing 
surface can be designed everywhere outside a small neighbourhood of the subsonic 
edges, L,T,, L,T,, where there is a vertical tangent. Alternatively the pressure 
distribution may be suitably modified near the edges so that the integral of 
expression (30) yields a solution for 0¢/0z which is regular at the wing tip edges 
(ie. the tip edges are then tapered). 


5. The Arrowhead Wing 


A wing of arrowhead plan form (bounded by L,OL,, T,TT, in Fig. 5) is 
considered when the edges are everywhere supersonic. The pressure on the wing 
surface is so distributed that the isobars are parallel to the leading edge and the 
pressure varies linearly with the normal distance from the leading edge. 


1 


Fig. 5. 


Let the angle of sweepback of the wing be - such that tan +<B. The isobars 
| are lines xcos7+ysin7=constant, the positive sign being taken for the left wing 
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tip. In a similar manner to the method indicated in Section 4 the final expression 
for the potential on the wing surface is 


o=C, (xcosr+ysinr’+D, (xcosr+ysinr), (49) 
or in characteristic co-ordinates, 


o= { C,/(4 B’) }{ «(Bcosz sinr)+ (B cos r+sinr) } 7+ 
+{ D,/(2 B)}{ 8 (Bcos7+sinr)}, (50) 


where C, and D, are known constants and when multiplied by cos7 correspond 
to C, and D, used in Section 4. 


09/0z at a point (x,y) in region L,OM,T, can be evaluated from expression 
(20) using the method of Section 4, or the problem treated as a rectangular wing 
whose leading edge is at an angle 4x—7; to the stream and result (45) used. 


Thus (0¢/0z).-+0= — /{ B? cos?z—sin? . (51) 


from which, using condition (5), the wing surface may be designed. 


0¢/0z at a point (2,8) inside the Mach cone from the apex is rather more 
difficult to evaluate and the method adopted will only be indicated. For a point 
in region OTM, the limits of integration of expression (20) can be divided (see 
Fig. 5): 


Oz/ 2-40 { (2—2’)(B— } 


3,+3,+5, 


(52) 


{ (2—2’)(B—B’)} 


=,+2, 


from which the final expression is found as 
(8¢/0z):- +0= — /{ B*? cos? r—sin? r }[2 C, (x cos sin r)+D,]+ 
+[(2.C, sin 2r)/=] (x? — B? y?}+ 
+{[4(2C,x cos 7+ D,) sin 7]/=}log[/(2By)/{/ (x— By)+ /(x+ By)}]+ 
+(2/{ B? cos? 7 — sin? r 2C; (x cos 7—ysinr)+ 
+D, } sin-'/{ (Bcos7+sin r)(x— By)/(x cos r—y sin r) } - 
— {2C,(xcosr+ysinr)+ 


+D, } sin~? { (B cos sin r) (x— By)/(x cos 7+y sin r)} 
(53) 
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As the x-axis is approached 09¢/0z becomes infinite, suggesting that pressure 
distributions with discontinuous gradients yield singular solutions. Since the 
singularity behaves like log y as y—>0, it is possible from equations (5), (51) and 
(53) to design the wing surface everywhere, except for a small neighbourhood of 
the line OT (see Fig. 5) along which the pressure gradient is discontinuous. A surface 
of reasonable shape, continuous with the theoretically determined surface, may now 
be chosen to complete the wing design across the neighbourhood of the singularity. 
This would give approximately the required pressure distribution outside the 
neighbourhood and a modified distribution inside. Alternatively by making the 
pressure gradient continuous across the x-axis in the theoretical considerations, the 
singularity would probably not arise. 


6. Further Problems 


Other pressure laws could be studied so that singularities which have arisen in 
the previous two cases do not occur. Another problem, which is of interest, comes 
from taking the isobars, obeying a similar pressure law as in the present case, to 
radiate from the vertices at the ends of the trailing edge of a triangular wing. This 
problem has been studied by the author and, although an answer could be obtained 
(by numerical methods, if necessary), it does not seem to admit a result in terms of 
elementary functions. 
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Longitudinal Waves in Beams 


R. E. D. BISHOP, B.Sc., Ph.D. 
(Ministry of Supply) 


Summary: The elementary engineering theory of longitudinal waves in beams 
allows for no dispersion, whereas Love’s corrected version of it predicts that very 
short waves are not propagated. A new correction is added which overcomes these 
defects. To do this, a simple theorem is proved concerning certain types of closed 
region. A method of improving Love’s formula and the present one is illustrated. 


1. Introduction 


When problems involving the flexural motion of beams are investigated, use 
may be made of the well-known simple theory (see, for example, Ref. 11, Section 
54). If greater accuracy is required (especially as regards short waves), two more 
exact versions are available. These are presented by Rayleigh’ and Timoshenko"? 
as corrected forms of the most elementary theory in which two of its known short- 
comings are rectified. As Timoshenko"® shows, his equation is remarkably 
accurate.* These three may be termed “engineering theories” and they owe their 
importance to the fact that the exact results of dynamical elasticity are far too 
complicated for any but the simplest problems. All of the engineering formulae 
are obtainable from elementary (strength of materials) reasoning and all tend to 
exactitude as the waves concerned become longer. 


The state of affairs is not quite as happy for longitudinal vibrations in beams. 
Here, again, exact theory is too intractable for most practical applications and it is 
necessary to turn to engineering theories. The simple theory is well known (see, 
for example, Ref. 11, Section 51) and it approaches agreement with the exact one 
for long waves. But its prediction of a fixed propagation speed (= /(E/ p)) for all 
wave lengths is contrary to the exact results and to what is actually observable, as 
found by Davies’. This simple theory has an obvious defect in that no regard is 
paid to lateral inertia forces, although lateral displacements must necessarily occur. 

A correction to overcome this is offered by Love? and his results show good 
agreement with those of the exact equations over a wider range of wave length than 
before. However, Love’s equation indicates that short waves are transmitted at 
vanishingly small speeds and this is neither compatible with the triie results nor 
with physical intuition. A less important point perhaps, but one worthy of 
consideration, is that this correction is obtained from Hamilton’s Principle and is 


*See also the work of Prescott™ and Christopherson), 
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thus unable to take its place with all the other equations in the average book on 
strength of materials. 


Notation 


P 
¢ (6,, 4) 


area of cross section 

area defined in Section 3 (see Fig. 2) 

radius of cylinder 

phase propagation speed 

phase propagation speed of long waves [= /(E/p)] 

speed of propagation of distortion waves in the infinite medium 
[= 

speed of propagation of surface waves 

Young’s modulus 


polar radius of gyration of cross section 
cylindrical co-ordinates (see Fig. 2) 


defined by equations (19) 

defined by equations (21) 

time 

displacements 

rectangular co-ordinates 

maximum values of y and z (Fig 1) 
normal stresses 

shear stresses 

wave length 

modulus of rigidity: one of Lamé’s constants [=E/{2(1+v)} ] 
Poisson’s ratio 

wave number (=27/A) 

density 


function defined in Section 3. 


2. New Theory 


The simplest theory of longitudinal motions of beams rests upon the 
assumption that, while vibrating, planes perpendicular to the axis remain plane. 
No restraint is allowed on the lateral motion of particles but the inertia forces 
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associated with it are neglected. Love’s correction is concerned with these forces 
and a correction for them is also devised here. 


Consider a prismatic beam whose cross-sectional area is A and whose centroidal 
axis is identified with the x-axis, and contemplate a thin slice of the beam of 


Fig. 1. 


thickness dx (see Fig. 1). Applying the law of motion to a small element dxdydz 
at P in the y and z directions, two of the equations of motion are obtained, i.e. in 
the usual notation, 


Ox dy 
Ores | OT ys 
dy Or’ 


N 
| 
> 
~ 


(1) 


The 
as a 


The 


as 


Sub: 
seco 


The 
as 


$0 


App 


or 


= 
| 
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These can be simplified if it is assumed that plane sections remain plane and that, 
as a first estimate, the displacements in the y and z directions are 


These are Love’s assumptions; but he omits consideration of shear stresses. With 
x as the modulus of rigidity, 


ov 
Try = ax vey ax 

ow (3) 
0 


Substituting into (1) and integrating the first equation with respect to y and the 
second with respect to z, 


2_ Ou = 2_y2 Ou 


The quantities y., Z. represent the maximum values of y and z corresponding to P 
as shown in Fig. 1; in general, they will be functions of the co-ordinates (y, z) of P. 


Now, at P 


ox (Tyy + 22) . . (5) 
so that 
+ [(ye Ze) (y +2" (Oe +2,7)—(y 
©) 
Applying the law of motion to the thin slice in the x direction then, 
2 
| dA | =Adx 
or? 
A 
OC 
A 
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On substituting from (6) into (7) and carrying out the indicated differentiation and 
integration, 


-(y? +Z *\| dA (8) 


A 


This equation may be taken as that of longitudinal motion of the beam, the functions 
ye and z, being introduced for any particular cross section. Obviously (8) would be 
a most difficult equation to work with for any but the simplest sections. However, 
it can be greatly simplified if use is made of the following mathematical result (which 
will be demonstrated in Section 3):— 


(y.? + =3 (y?+27)dA 


A A 


If this be accepted temporarily, the equation of motion may be written 
(10) 


where c, is the “ bar velocity” [= /(E/p)], k, is the polar radius of gyration of the 
cross section A and c, is the speed of distortion waves in an infinite elastic medium 


[= 


If, in the foregoing arguments, all stresses other than c,,, had been neglected, 
the elementary equation of long waves would have emerged. In point of fact, it is 
equation (10) with the two fourth derivatives omitted. 


Next, Love’s equation can be obtained by contemplating only c,.., 7,, and ¢, 
It is identical with (10) if the term involving 0*u/@x* is excluded, as may be seen b) 
tracing through the analysis. 


3. Mathematical Theorem 


Consider the area A’ which is bounded by the axes OY, OZ and a curve ( 
which is given by y=F (z), where F (z) is a single-valued function of z which decreases 
monotonically as z increases (so that F’(z)<0O for z>0). Suppose that C is 
represented by r=f(@) in polar co-ordinates where f(@) is single-valued in 6. Le 
P be some point in A’ with co-ordinates (r, 6) and let the perpendiculars from P to 
the axes intersect C as shown in Fig. 2, at P, and P, (with co-ordinates (r,, 4,) and 
(r,, 6.) respectively). Let the perpendicular to OY through P, intersect that to 027F 
through P, in P’ whose co-ordinates are (R, 9) (see Fig. 2). 
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Fig. 2. 
It is to be proved that, if dA’ is the element of area surrounding P, then 


R*dA’=3 | rdA’. 


A’ 
Since dA’=rdrd6 and the equation of C is r=f(@), it follows that 


6) 


A’ 0 0 0 
Now r,=f(@,) so that 
f (@,)cos 8,=rcos 6, 
and r,=f(@,) so that 
f(@,)sin 6,=rsiné. 


Therefore R? sin? 6, cos”6, + r? 


and | Reaa’= | 


A A A 
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That is, it remains to be shown that 


2/2 
r,2da’+ | r2da’=4 | rdA’= 
a’ a’ 


Consider the first of these integrals. Evidently, 


Now r=f (9,)cos 6,/cos @ so that, as r varies from 0 to f (8), 6, varies steadily from 
to 6. Thus, 


dr=[f' (6,)cos 6, —f (6,)sin 8,] 


so that 
2/2 6 
sec? 6 dé Uf cos 6, —f (6,)sin cos 6, dd, 
A’ 
/2 6, 


= | FOP F@,)sin 6, @,)c084,] cos 6, | sec? 


0 


That the last change of limits and of the order of integration is valid will be seen 
by reference to Fig. 3, the vertical axis being that of the function ¢ (6,, 6) where 


(9,,6)=[f [f @,) sin 4, — f’ (@,) cos 6,] cos 6, sec? 6. 


It will be noted that both double integrals represent the volume over the shaded area. 


Therefore, 
2/2 
r2dA’= (OF [f (6) sin 8 f (8) cos 6] sin 
A’ 0 


the suffix 1 being dropped. 
The integral 


/ 


| r2da’= 
0 


A’ 


1(6) 


do 
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(8,8) 


| 
| 
| 
from 
2 
ZZ 
seen 2 8, 
Fig. 3. 
can be treated in a similar way. Using the relation r=f(@,) sin @,/sin 6, this becomes 
rea. 
2/2 6 
| r,7dA’= cosec? 6 dé | Lf Lf (6,) sin 6, + f cos sin 6, 
A’ 0 0 
2/2 
= | IF @,)sin 6, +f (6,)c0s 6,] sin 4, cosec? dé 
0 65 


by the same sort of reasoning as before. Thus 


r,{dA’ = [f (9)]* Lf (6) sin 6 + f (6) cos 6] cos 6 dé, 
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the suffix 2 now being dispensed with. It is now clear that 


A’ A’ 0 
so that the theorem is proved. 


At this point a word might be said about the beam cross sections which are 
covered by this theory; for restrictions have been placed upon the curve C. The 
value of the corrections for lateral inertia which have been discussed seems 
questionable for sections other than those allowed in the proof. For example, no 
allowance can be made for the fact that the inner flange surfaces of an I-beam 
remain stress-free during longitudinal motion. However, it may happen that the 
empirical use of the equations is justified even for such sections as this. Of course, § 
no difficulty arises with square-, rectangular- or circular-cross sections. 


4. Comparison with Exact Theory 


The most convenient theory with which to compare these results is that of 
Pochhammer“? referring to a cylindrical shaft; but, even for such a simple section 
as this, the dispersion equation (governing wave length A and wave speed c) offers 
an unpleasant task of computation. However, Bancroft" gives a satisfactory table 


ed 1-0 B 


te) 


of nu 
wave 
the e1 
7 
where 
of the 
practi 
theort 
7 
by tr 
» 
\N These 
0-8 i I 
° 6 s/o the t! 
If thi 
D0 oo, 
0-3 INS secon 
° 2 4 6 4 10 12 14 16 18 “ 
as2A a-3A7° 
Fig. 4. 
288 


are 

The 
seems 
le, no 
‘beam 
at the 


ourse, 


lat of 
ction 
offers 

table 


LONGITUDINAL WAVES IN BEAMS 


of numerical values of c/c, and (effectively) a, a being the radius and 27/€ the 
wave length; this he does for several values of v. Taking v=0.3, his results are 
shown in Fig. 4 as curve A and this will be used as the standard of comparison for 
the engineering formulae. The curve has a horizontal asymptote corresponding to 
¢, the speed of surface waves. 


This exact theory gives rise to an infinite number of higher modes of vibration 
whereby nodal cylinders appear in the longitudinal motion. The analytical meaning 
of these has been examined by Cooper“? although they do not seem to have much 
practical importance (see Ref. 4). They have no counterpart in the engineering 
theories so that only the fundamental mode A need be considered. 


The dispersion equation appropriate to the simplified theories may be found 
by trying a solution of the form 


in the various equations of motion. This leads to the relations 


(elementary theory), ‘ . (12) 


(Love’s theory), ‘ . (13) 


(present theory). ; . (14) 


These are plotted as curves B, C and D in Fig. 4, respectively. For this purpose k,’ 
isreplaced by a?/2 and the value v=0.3 is used so that c,*/c,?=[2 (1+ »)]~*=0.385. 


Perhaps the most noticeable property of the curves is that D not only exhibits 
dispersion, but also has a horizontal asymptote of finite value. That of the exact 
theory corresponds to c/c,=c,/c,=0.575 while the one for equation (14) is 
c/¢,=c,/c,=0.621 when v=0.3, an error of 8 per cent. Simple theory is not likely 
to lead to any better agreement than this because surface waves find no place in 
the theory of strength of materials. 


One aspect of the new curve (D) may warrant further refinement; for while it 


isa great improvement upon B for all wave lengths, it is inferior to C for long waves. 


If this should be desirable (on account of the physical nature of the problem to be 


_handled by this theory), the method of Section 5 may be used. 


It is pointed out by Pochhammer'? that, if all higher powers of af than the 
second are rejected from the exact theory, the wave speed for long waves is given by 


(13a) 


c 
Co 
1 + 
Caley 
2 
Co” 
289 


R. E. D. BISHOP 


This coincides with Love’s theory, as may easily be found from equation (13).* [f 
this degree of approximation is assumed in (14), however, we find 


so that a new term appears (which is always less than half the one preceding it), 
5. Improvement of Corrected Theories 
In searching for a means of improving the engineering formulae, it seems logical 


to try to improve equations (2). They are obtained by integrating the relations 


ov 
=Ty,— V + 22) 


(15) 


0 
Ex 
under the assumption that o,, and o,, are small compared with c,, and can therefore 
be neglected. From this, the three stresses were calculated (equations (4) and (6)). 


Better agreement can be obtained by using the values of c,., v,, and ,, (found 
in (4) and (6)) in equations (15) so that improved expressions for the displacements 
v and w emerge on integration. If this is done, 


(16) 
—y*)zQ 
YY 3 ax +p (z,? —z yy ox’ 
oP 
Tw = 5 Ve 6y."y + 5y.*) fp or? 
4 — 2") (y.2-y Oe 
— Ou 2 p - 
Ze r 4 12 z — + 5Z6 ) or Ox? 


*It must not be concluded that Love’s theory is a second approximation of the exact one in 
any other sense, however. 
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where 


~2c,2 Loxdr? 


Proceeding as for equations (5), (6) and (7), the improved equation of motion is 


(19) 


12 \dxer? ox ox 
] 


(y* — 6y."y? + + — 62.727 + 5z.")dA , 


where R= 


| (21) 


Since no general result comparable to (9) can be found now, the particular case 
of a cylindrical shaft will be taken. For this cross section, 


Thus (20) becomes, here, 


[The newly added correction vanishes if the material is incompressible (i.e. if 
v=1/2)]. The corresponding dispersion equation is 
4 v2q? 


+ ae? + (2v? — | =0, 
(24) 


| and it is shown in Fig. 4 as curve E. 


There is no reason why this sort of improvement should not be made on Love’s 


| equation. Remembering that all shear stresses must now be disregarded, the 
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equations may easily be shown to be exactly as before except that all terms 
containing » and c,* are omitted. Thus, the general equation (comparable 


with (20)) is 


or? 0x? 0x? dr? 24c,? oxor* 
4 On 


and the dispersion equation is 


This last result is plotted as curve F in Fig. 4. 


OF axe 2 * 


(25) 


(26) 


(27) 


The curves show better agreement with the exact one. This is strikingly so 


with Love’s equation over quite a wide range of (long) wave lengths. 


This fact 


might be used if the distortion of a pressure pulse (as it is propagated) is being 


studied by consideration of the first few terms of its Fourier representation. 


6. Conclusions 


The engineering theories of longitudinal waves in beams have been examined 
in terms of the simple notions of strength of materials; a derivation of Love's 
correction has been given which is thought to be simpler than the original one (at 
least for the non-mathematically minded) and to give a clearer insight into the 
nature of the assumptions made. In doing this, a simple mathematical property of 
certain bounded regions has been demonstrated which might conceivably have other 


uses in the theory of second moments of areas. 


A new corrected form of the simple theory has been offered and its dispersion 
curve is better than those of the other engineering formulae for short waves. Where 
long waves are concerned, the new equation is inferior to Love’s version but is still 
better than the simple theory. The method of successive approximations has been 
used to improve the corrected formulae (both Love’s and the present one). 


There is a wide literature (both theoretical and experimental) on the problems 
of longitudinal impact. This is almost entirely based upon the simple theory and 
hence makes no allowance for the dispersion effects which undoubtedly occur and 
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which are suggested by dimensional analysis. Davies'‘) has used Love’s equation 
in treating such problems and he points out the anomalous results which spring from 
so doing on account of the zero phase velocity predicted by this theory for very 
short waves. The new theory should overcome this difficulty and open up a way 
of treating longitudinal impact with better accuracy than that offered by the simple 
theory. The exact theory is not a practical engineering tool. 
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A Note on Supersonic Wing Integral 


Equations in Unsteady Flow’ 


JOHN W. MILES 
(Fulbright Lecturer, Auckland University College, New Zealand. 
On leave from the University of California, Los Angeles) 


Summary: The integral equations for the pressure distribution on an oscillating 
thin wing having a prescribed velocity distribution in'‘a supersonic flow are developed. 


1. Introduction 


A plane wing, projecting on the region S of z=0, is placed in a supersonic flow 
of velocity U, directed along the positive x-axis, and executes a transverse motion 
of time dependence exp(iw?). If the co-ordinates are dimensionless, the linearised 
equation for the velocity potential is 


dyy + = + 2ikMo, — k?M*o, (x, y; Zz) . (1) 


(3 


where / is the characteristic length. 


2. Supersonic Trailing Edge 


Considering first the anti-symmetric problem where the trailing edge is every: 
where supersonic, the boundary conditions on the plane z=0+ may be cast in 
the form 


(x, y,0)=0, (x,y) not in (4 


*This work was done in 1948, at which time it was hoped to extend the very elegant methods 
available for the steady flow problem. All attempts to generalise these methods for wings 
of arbitrary plan form having proved nugatory, it is felt that the integral equations are 
worth recording. 


Paper received December 1951. 
[The Aeronautical Quarterly, Volume ITI, February 1952] 


294 


the 


wh 


ou! 


wh 
unc 
in 
the: 
car 
we 


As 
circ 
con 


Sl 

dis 

vir 

in 

pla 

|__| 


iting 


very: 
ist in 


(4a 


(4b 


ethods 
wings 
1s are 


SUPERSONIC WING INTEGRAL EQUATIONS IN UNSTEADY FLOW 


Equation (4b) is not actually valid in the wake, across which the potential is 
discontinuous, but suffices for the determination of the potential on the wing, in 
virtue of the (assumed) supersonic trailing edge. However, having the potential on 
the wing, the potential on the upper surface of the wake is given by 


¢ (x, x«>-x, (y) (5) 


where x, (y) specifies the trailing edge. 


In addition to the boundary conditions (4), it is implicit that ¢ must vanish 
outside the zone of action of S. 


A solution to equation (1) that reduces to a prescribed potential on z=0+, 
in consequence of Fourier’s theorem, and is anti-symmetric with respect to this 
plane is given by 


o(x,y,2=+ dy | dv | ae { ay x exp {i[u(x E)+y¥ (y—n)—A ] (E, n,0+) 
x 

(6) 

A= { v?-[(M— 1) n+ kM][(M +1) 4+ kM] (7) 


where < is a positive, real constant, placing the path of integration in the »-plane 
under all singularities therein, and the branch cuts in the v-plane are taken to infinity 
in the same half planes as the respective branch points from which they originate, 
these conditions ensuring .#(A)>0. Invoking the boundary conditions (4) and 
carrying out the » and »v integrations with the aid of a table of Fourier integrals, 
we obtain the integral equation 


fea (x,y) in (8) 
B? 
g(x, y)= =Ri [cos x>Bly|. . (9a) 
=0, x<Bly| . (9b) 


As written, the integral of equation (8) is improper, but this difficulty may be 
circumvented either by repeated integration by parts, or by adopting the well-known 
convention of Hadamard. 
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In the special case where the leading edges also are supersonic, 2 vanishes off 
the wing, and equation (8) may be inverted with the aid of the Faltung theorem to 
obtain the result first given by Garrick and Rubinow"’. For a rectangular wing 
of effective aspect ratio greater than unity the Wiener-Hopf technique is applicable”, 
In principle, an extension of this technique, as applied by Fox to the diffraction of 
sound by an infinite strip®’, renders possible closed form solutions for rectangular 
wings of arbitrary aspect ratio, but the resulting integrals are intractable, as is also 
the case if a solution in Mathieu functions is attempted after introducing co-ordinates 
of the elliptic cylinder. (The low aspect ratio rectangular wing may be so treated, 
however.) A workable but involved solution for the narrow delta wing may be 
obtained along the lines proposed by Stewartson™). Useful solutions for more 
general plan forms would appear to be feasible only with the aid of numerical 
techniques. 


3. Subsonic Trailing Edge 

In the more general case where the trailing edge of S is partially or totally 
subsonic the boundary condition (46) must be replaced by the requirement that the 
pressure, which is proportional to 


y (x, y)=9, (x, y,0+)+iko (x, y,0+) 


must vanish for z=0, (x,y) not in §. Equation (6) may be extended by replacing 
@ by —i(u+k)~*y in the integrand, and the resulting integral equation consequent 
on equation (4a) is 


‘8 
(k 
h(x y)= | exp| i x>Bly| . (14a) 
=0. x<Blyl . (14) 


The prospects of effecting a general solution to equation (13) for non-vanishing { 
are not encouraging. 
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The Theory of Aerofoils in Unsteady 
Motion 
J. R. M. RADOK 


(Department of Supply and Development, Aeronautical Research 
Laboratories, Melbourne) 


Summary: The theory of aerofoils in unsteady flow, which has made substantial 
progress in the last decade due largely to the ground work of Kiissner and his 
co-workers, is presented here in a form suitable for application in aeroelastic 
problems, particularly those concerned with the dynamic loads on aircraft arising 
from gusts. 

Exact expressions are given for the aerodynamic lift and moment for an 
oscillating aerofoil as well as for the case of arbitrary motion through disturbed air. 
The expressions for the latter case involve two special functions, generally referred 
to as Wagner and Kiissner functions. Exact values of these functions are tabulated 
together with useful approximations. The problem of a wing-tail combination is 
discussed and a method of solution indicated. The bibliography at the end of the 
paper lists the most important publications in this field. 


|. Introduction 

The loads induced in aircraft arising from dynamic causes, such as from 
encountering disturbed air, are becoming of increasing importance as design 
problems. An important factor in dealing with such problems is the behaviour of 
aerofoils in unsteady flow, a study which has made substantial progress in recent 
years due in the first place to the excellent theoretical ground work done by Kiissner 
and his co-workers about ten years ago. 


It is the purpose of the present paper to summarise the theory of the two- 
dimensional aerofoil in incompressible unsteady flow and to present it in a form 
suitable for direct application to aeroelastic problems associated with subsonic 
aeroplanes without sweepback. 


Several surveys already exist" * which deal with the earlier investigations 
into the aerodynamic forces acting on aerofoils in arbitrary motion. Hence the 
work done before 1936 has not been included here, particularly as the exact solution 
has only been obtained since then. The three-dimensional incompressible, as well 
as the two- and three-dimensional subsonic, problems have not yet found satisfactory 
solution. The supersonic problem in two dimensions®*~*") and in three dimensions 
for certain delta wings’? has been solved, but may be considered of minor 
importance for aircraft design at present. 


Paper received July 1950. 
[The Aeronautica! Quarterly, Volume III, February 1952] 
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The treatment of the present paper follows closely that of Kiissner:* who 
deduced his solution from that for the oscillating aerofoil. The solution of this 
latter problem will only be stated here since a number of independent deductions 
of the relevant results may be found in Refs. 7, 8, 37 and 38. 


The underlying assumptions are the same as those generally made in aerofoil 
theory, namely: 


(i) the aerofoil is assumed to be thin and to approximate to a flat plate, 


(ii) the constant forward velocity of the aerofoil is assumed to be large 
compared with the disturbances caused by its rotation, its lateral displace- 
ment and to the gust velocity, 


(iii) the disturbance is allowed for by adjustment of the boundary conditions 
on the aerofoil and the fluid is assumed to be irrotational everywhere except 
in the wake, 


(iv) the Kutta-Joukowski condition is assumed to ensure smooth flow at the 
trailing edge. 


Some of the results, given by Kiissner in 1940, had been obtained previously 
by himself and other investigators, but his final formula for the pressure distribution 
over the aerofoil and his method of deduction are different from all earlier work. 


Notation 
frequency parameter 
21 chord of aerofoil 
V constant forward velocity of aerofoil 
w=ivl/V reduced frequency 


x  chordwise co-ordinate of a point on an aerofoil such that x= ¥! 
corresponds to leading and trailing edges respectively (see Fig. |) 


W _ vertical velocity 
6 cos~'(—x) 
t time 
variables of integration 
II pressure difference 
P,,@, constants in equations (1) and (2) respectively 
Hi, (2, Hankel functions 
C Theodorsen function 
T 2C-1 
p density 


wi 
W 
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lateral vertical displacement of aerofoil 
angular displacement of aerofoil 
distance of point of rotation from centre of aerofoil 


distance of point about which moment is taken from centre of 
aerofoil 


lift 

moment 

instantaneous pressure difference 
transient pressure difference 

Vt/l, non-dimensional time parameter 
vertical velocity of points on the aerofoil 
Wagner’s function 

gust velocity 

amplitude of W 

Bessel functions 

Kiissner’s function 


values of 6 between 0 and = (see equation (21)) 


special function describing growth of instantaneous pressure 
difference due to gust 


bound vorticity 


distance of centre of tail section from centre of wing section 
tail chord 
U,+1 
C1 distance between the centre of the wing section and a point in the 
wake. 


2. The Harmonic Oscillating Aerofoil 


Familiarity will be assumed with the concepts of bound and free vorticity and 
with the deduction of some of the basic results, referring to the oscillating aerofoil, 
which follow. 

Let »v be the frequency parameter of the motion 

w=ivl/V, the reduced frequency 
V, the constant horizontal, i.e. forward, velocity 
W, the small vertical velocity 
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x=-—cos 4 the co-ordinate of points on the aerofoil of chord length 2. 
so that 6=0, 7, ie. x= ¥1 correspond to the leading and 
.iling edge respectively (Fig. 1). 


+ | 
Ww 
Fig. 1. 
Write W (6, t)=Ve™ [P, +22 cos n6] (1) 
1 


where the P,, are constants. The corresponding pressure difference distribution can 
be assumed in the form 


Il (6, t)=p V? e*[2a, cot 46 +4 Xa, (1/n)sin n9] (2 
1 


where the a, have to be determined in terms of the P,. Using the condition of 
tangential flow on the aerofoil, the following general relations are found between 
the P,, and a,: 

a,=4(1+T)(P,—P,)+P, 


nP,,- n>1 


(3) 
where H, |.) (—iw), H, (.) (—iw) are Hankel functions. It is easily verified that 
T=2C-1, 
where C=H,,.)(—iw)/ { Hy) (—io) +i Ay.) (—iw)}. . (4) 


is Theodorsen’s function”. The relations (3) were first stated by Kiissner'”’ in 1936; 
their detailed deduction was given by K. Jaeckel"''”) in 1939. 


Thus it is seen from (3) that with a harmonic oscillating aerofoil, three distinct 
stages lead to the determination of the pressure difference distribution (2) over 
the aerofoil; these are: 


(a) Calculation of the P, as Fourier coefficients of the given W (8, t) 
(b) computation of the a,, using (3) 


(c) summation of the series (2) after substitution for the a, in terms of the P,. 
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Kiissner succeeded in performing some of these steps once and for all and in 
deducing a formula for I1 (6,1) which involves only a single integration over the 
aerofoil after substitution of the known W (6,1). The indirect mathematical 
reasoning, leading to this formula, is related to the idea of the Poisson integral; 
it may be found in Ref. 5. The final result, determining the difference in pressure 
between points above and below the aerofoil for harmonic oscillations, is given by 


pV 2 sin 6 1 —cos (6+ #) 
II (, t) | - + sin? log of 
= J Lcosé—cos# -1—cos 
0 


+ { 1+cos#+(1—cos 9) T (—iw) } cot4@ | W (I, t)dd . (5) 

Next the lift on an oscillating aerofoil, as given by Theodorsen'’? [cf. equation 

XVIII, p. 419], will be obtained, in order to demonstrate the application of (5). 
For this purpose let 


W (8, t)- { h—(a+cos 9) 1z+V2 } 
—{ h, (a+ cos #) 1a, } iv + Vx.) (6) 
where A,, z, are the amplitudes of the lateral and angular displacements fA and z, 


2/ is the chord of the aerofoil and a is the distance of the rotational centre behind the 
centre of the aerofoil (Fig. 2). 


Fig. 2. 


Substituting (6) in (5) and integrating with regard to 3 we obtain 
(6, — pV [2sin@ { 21z+(l/V)h—al? 4} + 
(7) 
where, as before, dots indicate differentiation with respect to the time f. 


Using (7), the lift L acting on the aerofoil at time ¢ is given by 


I 


11 (6, t)dx=I 


II (6, t) sin 6 dé ; : ; (8) 


o—, 
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ie. L=— pV h—-al a} + 3] 


and by (4) 

This is Theodorsen’s result after putting V=v and /=b and remembering that in 
the present case the flap angle 8=0, i.e. that there is no aileron. 

Similarly the moment about the point b is given by 


0 t 
i.e. for b=a: 


M= — [x (4—a) V124+ (4+?) at lh] + 
+2p VP x(a+4)C { Vat+h+1(4-a)2 } (ll) 


(cf. Theodorsen’) equation XX, p. 419) 


3. The Arbitrarily Moving Aerofoil 


An inspection of (5) shows that the expression for II (6, t) can be split up into 
the two parts 


2 sin 


II, (6, +cos cot 4 6 + 
cos #—cos # 


0 


+ 


ol en sinddd (12 


and IL, (6,1)= (pV T (—iw) cot $4 | (1 cos 9) W (0, 


0 


where in (12) use has been made of the relation 


which holds for the harmonic oscillating aerofoil. 


Kiissner calls II, and II, the “instantaneous” and the “transient” pressure 
difference. Remembering that T(—iw) originates from the presence of the wake, 
i.e. from the history of the motion, the suitability of these terms becomes obvious. 
Further it is clear that (12) must also hold for arbitrary motion, while (13) is only 
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THE THEORY OF AEROFOILS 


IN UNSTEADY MOTION 


valid for harmonic oscillations. The remainder of this section will be devoted to 
the task of deducing an expression for II,, applicable to an arbitrarily moving 
aerofoil. 


Kiissner achieves this object simply by writing down a general formula for IJ, : 


) W (i.e 
cot | (I cos #) | Us (s a9, (s>0) (15) 


Co 


IT, (6, s) 
0 0 
where W, is the vertical velocity component of points on the aerofoil, s=Vt// is a 


dimensionless co-ordinate along the flight path (Fig. 3) and U, is a function which 
is to play a similar role to that played by T and which has still to be determined. 


START OF 
AEROFOIL WAKE DISTURBED MOTION 
sel $ 3 2 
Vt 
Fig. 3. 


It is found by using the fact that (15) must obviously also hold true for a harmonic 
oscillating aerofoil, in which case it must be equal to (13) provided that s—>oo. 
The latter condition is necessary since the function T is obtained on the assumption 
of an infinite wake. 


0 
or by inversion of the transform, 
U, (s)= 5= |T do, s>0 . (16d) 
| 


where the path of integration passes the origin on a small semi-circle on the positive 
side of the real axis. 


In this case, the function U, (s) is to take account of the fact that any point 
4 of the aerofoil has the vertical velocity component W, (6,5) measured with respect 
to a zero component in steady motion. In a similar manner, by considering partial 
motion in the fluid and neglecting the effects of any vorticity in the fluid due solely 
to this motion, the transient pressure difference is obtained, caused by a “ vertical 
gust velocity” W, (s) acting in the positive (downwards) direction. Here 


s+l 


0 
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where the upper limit (s+ 1) takes account of the fact that the leading edge enters 
the “gust” ahead of the centre of the aerofoil, from which the s co-ordinate is 
measured (Fig. 3). The function U., is again determined by reduction to the 
oscillating case, which now involves substitution of 


W (3, t)=W, exp (os + cos W, exp (ivt — ox) 
in (13) and of (s)=W, 


in (17), i.e. a stationary field and an oscillating aerofoil and a harmonic oscillating 
field are considered respectively. Equating the resulting expressions leads to 


s+l1 


U, T (—iw)e*=0 | U, do for s> (18a) 
0 


or letting s—>oo and inverting the transform 


[J, (—iw)—iJ, (—iwv)] T (—iw)e dwo/wfor s>—1 . (18d) 


— 18 


U, (s)= 4 


Tl 


where J,, J, are Bessel functions and the path of integration passes the origin on a 
small semi-circle on the positive side of the real axis as in (165). 


By the help of the functions U, and U, it is now possible to write down the 
general expression for the pressure difference at any point @ of the aerofoil at time 
t=si/V : 


II (6, S) x { (1 + COS cot 3 6 cos 6 — cos 


0 


1—cos(@+7) 


s+1 
} Ow, 7 
+(1—cos 9) cot 4 {u, (s— U,(s—c) do \d9. 
0 
(19) 
Here use has been made of the identity 


0 


and of the fact that U,(s)=0 for s<0. The latter statement is obvious, since it 
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has been assumed that disturbances ahead of the leading edge have no influence on 
the motion of the aerofoil. By the same reasoning U,(s)=0 for s<~— 1. 


This expression (19) is Kiissner’s final result for the pressure difference 
distribution over an aerofoil in arbitrary motion. The expression involves a certain 
amount of vagueness since in actual fact W, also varies over the aerofoil, i.e. it is 
a function of s and 6. This will only affect the instantaneous part, since the transient 
part already involves an integration over the aerofoil and its wake. This difficulty 
can be overcome by a special study of the instantaneous pressure difference (12) 
due to a unit change in W,, travelling across the aerofoil. Thus 


for 0<0<9, let W,—W,, 
ct 
cos@—cosv’ ar cos @—cos 
and for Y<id<r, let W,=0, 
Substituting (21) in (12) and letting ¢—>v we obtain 
lim IL, (6, 1)= pV Jcot pV (22) 


where s= —cos®@ is the co-ordinate of the sharp gust-front on the aerofoil and 


P(s) sin ‘st (1—s*} for | 


(s)=1 for s>1 


It is easily seen that in the limit as ¢—>v the gust velocity W, given by (21) 
corresponds to a step change W, at the point —cos@ and that hence, by (72), the 
function P (s) describes the building up of instantaneous pressure as a unit nge 
in the gust velocity traverses the aerofoil. 


In the case of an arbitrary gust W, (s), the instantaneous pressure difference due 
to it can now be expressed, using the function ®, as 


s+1 


II, (6,s)=pV cot (s—c) 


0 
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Using (24) and (20) the formula (19) can now be rewritten 


II (6, s)= {a +cos #)cot4 6+ 
0 


2 sin 6 1—cos(@+¥) 
cos §—cos 1—cos(@—9)° as 


+(1—cos #) cot 46 {u, (s—c) r) + ky (s—0) do | as 


0 


} W, (3,5) + 


(25) 


Finally note that in (25) the integral with respect to « may also take account of 
discontinuous changes in the velocity functions. 


4. Arbitrary Motion of a Wing-tail Combination 


where y (6,5) is the bound vorticity distribution over the aerofoil, it is seen from 
(25) that in the case of arbitrary motion 


{a +cos #)cot4 6+ 


0 


2 sin cos(#+#) 
+sin?# log 
cos 6—cos # -cos (@—#)* as 


s+1 


(3, s)+ 


(28) 


Thus the bound vorticity distribution for any particular motion of the aerofoil is 
easily found. 


When considering a wing-tail combination, in addition to the forces due to 
the motion of the tail section and the gust, there will be effects resulting from the 
presence of the wing. The former forces will be given again by an expression of the 
form (25), the only difference being that (s—) must be replaced by (s—«)/£, where 
2é1 is the chord of the tail-section. Further account must be taken of the fact that 
the gust velocity function W, (c) will not act before o attains the value A—£ where 
Al is equal to the distance between the centres of the wing and the tail sections 
(Fig. 4). 
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The general problem for the effects of the wing and its wake on the aerodynamic 
forces on the tail is very much more involved. It is customary to assume that the 
chords of the wing and tail sections are parts of the same straight line which contains 
their wakes. Although this assumption is not essential it simplifies considerably the 
expression to be found later for the velocity induced by the vortex system of the 
wing on the tailplane. Once this induced velocity has been determined it can be 
added to the gust velocity W,(c) in the expression discussed in the preceding 
paragraph. Another assumption, which is not in agreement with actual facts, 


WING re TAIL 


| r END OF WAKE 
OF WING 


2 | O 


Fig. 4. 


concerns the effective incidences of the wing and tail, which have been assumed 
equal. This should not influence the final result to any great extent, since only 
changes from the equilibrium condition, i.e. from the condition of steady motion, 
are important. 


In order to determine the induced velocity on the tail it is necessary to determine 
the free vorticity distribution over the wing section and its wake. It can be found 
from the bound vorticity in the following manner, since 


Oy Oyds V dy (6,5) 
is the variation of the bound vorticity y at time t=/s/V at the point (—/ cos 4). 
The free vorticity at any point (—/cos4@,) on the aerofoil may be defined as being 
equal in magnitude and opposite in sign to the accumulation of the bound vorticity 
variations generated between the point under consideration and the leading edge. 
In this connection attention must be given to the fact that these variations move along 
the aerofoil with velocity V so that the variations constituting the free vorticity at 
(—Icos@,) at time Is/V originated at the earlier time (//V { s—(cos@—cos@,) } at 
the point (—/cos 4). Remembering that 6=0 corresponds to the leading edge, we 
have for the free vorticity 


(cos 6—cos 6,)] (—/ sin 6) dé 


6 
0 
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In particular, the free vorticity at the trailing edge is given by 


5, 6,5 cos 6— 1)sin 6 dé 
0 


and hence at any point (¢/ of the wake by: 


07 


0 
where (/ is the distance between the centre of the wing section and the point in the 
wake under consideration. 


Using these expressions for the free vorticity the velocity induced by the wing 
vortex system on the tail can be written as 


V Oy 
cos 6—()sin 6 dé — 


0 
1 
6 
4 sin6,d6, dy l y(6,s) . 
| — | —[6,s—(cos 6—cos sin 6 dé- sin 6 dé 
2x} cos 6,-¢ J as | J] 2x J cos 
0 0 0 


(33) 


where A—-€XCSA+€ is the range of the tail section. 


As already mentioned, this induced velocity can be added to the gust velocity 
term in the expression (25) after modification for the tail section. 


5. Remarks on the Special Functions 


It is seen from the final expression (25) for the pressure difference distribution 
that the functions U,, U, and k, are as fundamental to the problem of an aerofoil 
in arbitrary motion as is the function T (see (3)) to the problem of an oscillating 
aerofoil. These functions have been investigated by different authors who also spent 
considerable time on their tabulation. The formulae (16b) and (185) had been 
deduced previously by Kiissner) and I. E. Garrick"*'*). In their present form they 
are of little value for the determination of the actual values of the functions. 
Detailed analytical investigations of U,, U, and k, can be found in Ref. 6 (Section 
VI) and in a paper by L. Schwarz"®. 


In analogy with (26), i.e. 
k, (s)=U, (s)+® (s) s>-1 


308 


rel 


it 
V 

TI 

mi 
an 
ar 
ap 
ca 
wl 
fal 
R. 

|_| 


32) 


the 


ing 


33) 


ity 


34) 


THE THEORY OF 


AEROFOILS IN UNSTEADY MOTiON 


it is convenient to introduce a function k, (s) defined by 


k, (s)=U, (s)+1 for s>0 (35) 
The analogy will be still clearer when reference is made to (23), for it is seen then that 
k, (s)=U, (s)+1, k, (s)=U, (s)+ 1 for s>1. . (36) 


The function k, (s) was first calculated by H. Wagner"*’. Kiissner’ deduced 
the function k, and tabulated it together with k, in 1936. However, he made a 
mistake in the determination of k, which was first pointed out by von Karman 
and Sears"” in 1938. In 1940 Kiissner gave corrected values of k, and k, which 
are reproduced in Tables I and II towards the end of this paper. While these tables 
are sufficient for use with any approximate integration formula, there are also 
approximate expressions for U,, U, and k, which have been given by various authors. 
Some of these will now be given here and conclusions with regard to their accuracy 
can be drawn by comparison of the relevant columns of the Tables I and II in 
which some of these expressions have been tabulated. 


In chronological order the approximations are: 


4 


I. E. Garrick (1938, Ref. 13): k, (s)=2 a 


(37) 
Garrick’s approximation to k,(s) was incorrect due to the use he made of 
faulty results drawn from Kiissner’s paper’. 


I. E. Garrick (1939, Ref. 14): 


k, (s)= + (38) 


| 


This expression is the result of substituting (37) into the following formula, 
relating U, (s) and U, (s): 


+1 
1/2 
U, (s)= (s—x)(4**) ax. ‘ . (39) 

R. T. Jones (1939, Ref. 18): 
k, (s)=2 —0.33 994555 — 0.67 e- 9-35 : . (40) 

R. T. Jones (1940, Ref. 19): 
k, (s)=2 — 0.472 9-958 G+) 1,026 9364 6+) 0.342 e-24264D (41) 
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W. R. Sears and B. O. Sparks (1941, Ref. 30): 
k, (s)=2- 1.0 29-130 G+) _ 10e (42) 


W. Prichard Jones (1945, Ref. 20): 
k, (s)=2 —0.33 9-415 — 0.67 9-525, (43) 


In Ref. 20 the author also confirms Kiissner’s computations of k,(s) and 
reprints his table. It should also be noted here that some authors use 7k, (s) 
and <k, (s) instead of k, (s) and k, (s); while W. P. Jones denotes the new functions 
again by k, (s) and k, (s), R. T. Jones uses the symbols C;.(s) and Cr, (s). Further, 
in some places the argument of the functions U, (s) or k, (s) starts from s=0 instead 
of s=—-1. 


Finally certain physical interpretations of the functions U, (s) and U, (s) will 
be given, which will be of particular value in connection with the various 
applications of the theory in engineering problems, since that is the way in which 
some authors introduce these functions. Thus U, (s) describes the growth of the 
transient pressure difference following a sudden unit change in the vertical velocity 
component of the aerofoil itself. The function U,(s) takes account of the growth 
of the transient pressure difference subsequent to the entrance of an aerofoil into 
a sharp edged unit gust, while the function k,(s) represents the growth of the 
total pressure difference following a sudden unit change in the gust velocity. (For 
® (s) cf. remarks following (23).] On the basis of these interpretations the formulae 
(15) and (17) can be deduced by some physical reasoning similar to that used before 
for obtaining the expression (24) for the instantaneous pressure difference subsequent 
to a sharp gust. 


6. Expressions for Aerodynamic Forces Used in Applications 
Note: In the whole of this Section’ denotes differentiation with respect to s. 


The final formula (25) gives only the pressure difference distribution over an 
aerofoil in arbitrary motion. One further integration over the aerofoil is necessary, 
therefore, in order to obtain expressions for the lift and aerodynamic moment acting 
on the aerofoil [see (8) and (10)]. Thus, assuming that the aerofoil is rigid and is 
subject to vertical and angular displacements h and 2, where the latter occurs about 
an axis inside the aerofoil, and to changes in the vertical velocity component of the 
fluid, we can write for W, (3,5): 


W, (8,s)= +Va—(a+cos 91S | — [ b+ V2—(a+cos 


(44) 
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(see Figs. 2, 3). Substituting (44) and an arbitrary W,(s) in (25) and integrating 
with regard to # gives 


II (6,s)=pV [cot46 { A-4B+C)} +2sin6{ A’+B)} —sin@cos6B’] (45) 


where A(s)= h’+laz—ala’}, B(s) V 2’ 
| 
C(s)= U, (s—o) { A'+4.B’ } +k, (s—o)W,/]do. 


0 


‘ . (46) 


Substituting (45) in (8) and (10) and integrating with respect to 6 over the aerofoil, 
the exact expressions are obtained finally for the variations in lift L and moment M 
from the equilibrium position of steady motion. 


s+l1 
ptr V2 +h” a} —ala” + [ o){h” +12’ +12” (4—a)}do]+ 


0 


pr lV } k, (s—o)W,' do . (47) 


0 


dM=prPV[— {b+}}{ A-4B+C} —b{ +B’/8] 


+1 
pr IV2[— {b+4}{ U,(s—o){h” +12’ +(4—a) la’ }do } 
0 
s+1 


{ ) pr PV (b+4} [ o)Wide (48) 


0 
or using (35), i.e. 


U, (s)=k, (s)—1, s>0, 
s+1 
dL= —prV? [ +I2’—ala” + k, (s—o) { h” +12’ +12” (4-a)} deo'| 


0 


s+ 
+prlV | k,(s—c)W,'do. . (49) 
0 
dM = px (b/I)h” (4—b)+2” (4 +ab)- 
s+1 
~(6+4) | ky + +2" do | 
0 
s+1 
=pt PV {b+4} [ k.( o)W,' do ‘ (50) 
0 
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Note that as before al is the centre of pitch of the aerofoil and b/ is the point about 
which the moment is taken. 


In the case of a wing-tailplane combination similar expressions will be obtained, 
the main difference arising out of the term (33), which is to be added to W, (s/¢), 
the time lag, due to the position of the tailplane behind the wing section, and the 
modification to be made to W, (9,s/£) as given by (44). In actual fact W, (#,s/£) 
becomes in this case 

«8 V 
{ h’+la+(A—a—Ecos 9) Iz } 


« 


where, as in Section 4, the chord of the tailplane is given by 2¢€I. 

In many applications the co-ordinate s along the path of the aerofoil is measured 
from the leading edge instead of from the centre of the aerofoil so that the upper 
limits of the integrals entering into the expressions for the lift and moment become 
s instead of (s+ 1). 

It is often more convenient to present the formulae (49) and (50) in a different 
manner, the principle of which is to divide the lift force into localised components 
and to state with each of these their point of application so that the total lift and 
the moment about any point b/ of the aerofoil can be written down. By this method 
of representation the formulae given previously take the form (positive forces acting 
downwards): 


(a) Inertia forces 
—ptV? { h”—alz” } at 4-chord point 
— pt lV? a’ at }-chord point 


(b) Wagner force 


— pr lV? | k, (s—o){ at }-chord point 


0 


(52) 
(c) Kiissner force 
| k,(s—o)W,/ do at }-chord point 
0 


(d) Pure couple 
— pr PV? 2” /8. 


This classification of the various lift components is self-explanatory, the first 
heading being due to H. Glauert'*’), who gave it a physical interpretation. 

Naturally the representation (52) is not unique, but the one given here is that 
mostly occurring in applications. Further it was pointed out at the end of Section 
4 that the expressions for the Wagner- and Kiissner-forces are such as to take 
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account of discontinuous disturbances. In applications, the discontinuities often 
occur at the beginning of the disturbances (sharp-edged gusts) for which reason the 
terms just given for the Wagner- and Kiissner-forces are augmented respectively by 


— pt lV? k, (s){ A’ k,(s)W, (0) 


where h’ (0), 2(0), 2’(0) and W, (0) refer to initial discontinuous changes from the 
undisturbed state. 


In References 21, 24 and 28 the aerodynamic forces have been presented in a 
similar manner as in (52). It should be noted that in the first of these references the 
centre of rotation is the quarter-chord point (a= -—4). In the second reference'*’ 
the centre of rotation is at the three-quarter chord point (a=4). In the last paper'?*® 
the inertia term due to the angular acceleration is absent. 


Finally, mention should be made of two papers by von Karman and Sears"” 
and by Sears alone because they have given an independent deduction of Kiissner’s 
results and pointed out a mistake made by Kiissner® in 1936. These authors, 
defining aerodynamic lift and moment, respectively, as the rate of change of the total 
momentum and of the total moment of momentum of the aerofoil and wake, based 
their treatment on circulation theory which Kiissner avoided intentionally in his 
last paper’. The second paper by Sears is of special interest since it opens up new 
methods for the study of the special functions U, and U,. 


7. Three-dimensional Effects 


For the sake of completeness a few remarks will now be made about three- 
dimensional effects, taken into account by various writers. The general literature on 
three-dimensional aerofoil theory is larger than that on two-dimensional theory, but 
as mentioned in the introduction, it is still far from complete. For this reason all 
investigators of gust load problems have restricted themselves to the use of the 
results of the two-dimensional theory, but it is customary to make allowance for 
finite wings by replacing the factor 2 in the expressions for the aerodynamic forces 
and moments by the actual slope of the relevant lift coefficient curve in steady flight. 
However, this procedure may lead to errors with wings of low aspect ratio. There 
is also another factor, mostly taken into account by consideration of three- 
dimensional effects, which is the downwash behind the wing and its effects on the 
tailplane. The procedure is similar to the one adopted in Section 4, the main 
difference being that account is taken of an inherent instability in the vortex sheet 
forming the wake. 


Among recent investigators of three-dimensional aerofoil theory the most 
important are H. G. Kiissner“", L. Schwarz”, W. Prichard. Jones?” and E. 
Reissner“*-**), most of whom consider only the case of a harmonic oscillating wing. 
It should again be possible to deduce the arbitrary case from the oscillating one, 
once a satisfactory solution has been found for the other. 
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TABLE I 
WAGNER’S FUNCTION 


| | | Approximations to k, (s) 


U, (s) k, (s) | Equation | Equation Equation 


0.024 39 | 1.02439 | 1.02440 | 1.02130 | 1.022 45 
0.2 | 0.04764 | 1.04764 | 1.04762 | 1.04201 | 1.044 24 
0.3 | 0.06982 | 1.06982 | 1.06976 | 1.062 14 | 1.065 36 
0.4 | 0.09103 | 1.09103 | 1.09090 | 1.081 71 | 1.085 87 
0.5 | O11133 | 1.11133 | 12 | 110074 | 1.105 77 
06 | 0.13078 | 1.13078 | 1.13044 | 1.11926 | 1.125 06 
0.14945 | 1.14945 | 1.14894 | 1.13726 | 1.143 80 
0.167 38 | 1.16738 | 1.16666 | 1.15476 | 1.161 98 
09 | 0.18462 | 1.18462 | 1.18368 | 1.17178 | 1.179 62 
10 | 0.20121 | 1.20121 | 1.20000 | 1.18833 | 1.196 74 
11 | 0.21719 | 1.21719 | 1.21568 | 1.20443 | 1.213 35 


| Ref 6) | Re. 6 | G7 | | 43) 
0.0 | 0.00000 | 1.00000 | 1.00000 | 1.000 00 | 1.000 00 


12 | 0.232 60 | 1.23260 | 1.23078 | 1.22009 | 1.229 48 
13 | 0.24747 | 1.24747 | 1.24528 | 1.235 33 | 1.245 14 
14 | 0.26182 | 1.26182 | 1.25926 | 1.25014 | 1.260 35 
15 | 027569 | 127569 | 127272 | 1.26456 | 1.275 10 
1.6 | 0.289 10 | 1.28910 | 1.28572 | 1.27859 | 1.289 42 


1.7 | 0.30207 | 1.30207 | 1.298 24 | 1.292 23 | 1.303 34 
18 | 0.31463 | 1.31463 | 1.31034 | 1.305 51 | 1.316 85 
19 | 0.32679 | 1.32679 | 1.32204 | 1.318 42 | 1.329 96 
2.0 | 0.33858 | 1.33858 | 1.33334 | 1.33100 | 1.342 70 
2.1 | 0.35001 | 1.35001 | 1.34426 | 1.34323 | 1.355 06 
22 | 0.36109 | 1.36109 | 1.354 84 | 1.355 14 | 1.367 08 
2.3 | 0.37185 | 1.37185 | 1.36508 | 1.36673 | 1.378 75 
0.382 30 | 1.38230 | 1.37500 | 1.37802 | 139008 | 
0.392 45 | 1.392 45 | 1.384 62 | 1.388 99 | 1.40110 | > 


N 
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TABLE I (continued) 


IN UNSTEADY MOTION 


U, (s) 
(Ref. 6) 


0.4391 
0.4799 
0.5159 
0.5479 
0.5764 
0.6020 
0.6251 
0.6460 
0.6650 
0.6824 
0.6983 
0.7128 
0.7263 
0.7386 
0.7501 
0.7706 
0.7883 
0.8296 
0.8733 
0.9003 
0.9183 
0.9405 
0.9535 
0.978 12 
0.995 903 
0.997 973 
1.000 000 


| 
| 
| 


k, (s) 
(Ref. 6) 
1.4391 
1.4799 
1.5159 
1.5479 
1.5764 
1.6020 
1.6251 
1.6460 
1.6650 
1.6824 
1.6983 
1.7128 
1.7263 
1.7386 
1.7501 
1.7706 
1.7883 
1.8296 
1.8733 
1.9003 
1.9183 
1.9405 
1.9535 
1.978 12 
1.995 903 


1.997 973 | 


2.000 000 | 


Equation 
(37) 


1.428 58 
1.466 66 
1.500 00 
1.529 42 
1.555 56 
1.578 94 
1.600 00 
1.619 04 
1.636 36 
1.652 18 
1.666 66 
1.680 00 
1.692 30 
1.703 70 
1.714 28 
1.733 34 
1.750 00 
1.789 48 
1.833 34 
1.862 06 
1.882 36 
1.909 10 
1.925 92 
1.961 538 
1.992 064 


1.996 016 | 
2.000 000 | 


Approximations to k, (s) 


Equation 
(40) 


1.439 70 
1.484 12 
1.523 11 
1.557 41 
1.587 65 
1.614 39 
1.638 09 
1.659 17 
1.677 96 
1.694 79 
1.709 91 
1.723 53 
1.735 86 
1.747 06 
1.757 27 
1.755 23 
1.790 53 
1.825 79 
1.865 51 
1.893 83 
1.915 64 
1.946 53 
1.966 08 


1.996 513 | 
2.000 000 | 


2.000 000 
2.000 000 


Equation 
(43) 


1.451 65 
1.495 50 
1.533 63 
1.566 86 
1.595 90 
1.621 35 
1.643 74 
1.663 50 
1.681 00 
1.696 58 
1.710 49 
1.722 96 
1.734 22 
1.744 41 
1.753 69 
1.769 96 
1.783 84 
1.816 08 
1.853 55 
1.881 38 
1.903 50 
1.935 99 
1.957 52 
1.994 531 
2.000 000 
2.000 000 


| 2.000 000 
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| 35 | 
| | 40 | | | 
| 45 | | 
5.0 
6.0 | 
| | 6.5 | | | 
| | 7.0 
| 
| 85 | 
| 60 | 
9.5 
10 
12 
15 
20 
25 
30 
40 
50 
100 
|_| 
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KiiSSNER’S FUNCTION 


+0.1 


U, (s) 
(Ref. 6) 


0.000 00 


0.000 05 
| 0.000 27 
| 0.000 73 
| 0.001 51 


0.002 63 
0.004 15 
0.006 12 
0.008 56 


0.011 53 
0.015 08 


0.019 25 


| (Ref. 


k, (s) 


0.000 
| 0.282 


0.396 
0.481 


| 0.551 


0.611 


| 0.664 


0.712 
0.756 


| 0.796 


0.024 10 | 


0.029 70 


| 0.036 11 
| 0.043 45 
0.051 83 | 
0.061 42 | 
0.072 47 | 


0.085 44 
0.101 62 
0.194 2 
0.270 3 
0.334 4 

l 


| 0.389 


0.833 
0.867 
0.900 
0.930 
0.958 
0.985 
1.011 
1.035 
1.058 
1.080 
1.101 
1.194 
1.270 
1.334 
1.389 


6) 


00 


| 0.000 00 
| 0.282 31 
| 0.395 82 
| 0.480 51 
| 0.549 81 


0.609 00 
0.660 75 


| 0.706 66 
| 0.747 79 
| 0.784 83 


0.818 31 


| 0.848 60 


0.875 97 


| 0.900 63 


0.922 73 


| 0.942 33 


0.959 48 
0.974 13 


| Equation 
(38) 


Equation 
(41) 


0.000 00 | 0.160 00 
| 0.272 
| 0.368 
0.450 


0.282 
0.396 
0.481 


0.551 


0.611 
0.664 
0.712 


| 0.756 


0.796 


| 0.833 


0.986 15 | 


0.995 18 


(s)= 1 
for 


0.867 
0.899 
0.930 
0.958 
0.985 
1.011 
1.035 
1.058 


| 1.080 
1.000 00 | 


1.100 
1.192 
1.267 
1.329 
1.381 


35 
10 
27 
33 
66 
89 
73 
32 
36 


78 
99 
24 


0.521 


| 0.584 
| 0.639 


0.688 


| 0.733 
| 0.773 
34 


0.811 
0.845 
0.878 


0.908 
0.936 
0.963 
0.989 


25 | 


09 | 


98 
60 
12 
00 


20 


1.014 
1.037 
1.060 
1.081 
1.177 
1.258 
1.327 
1.386 


92 
70 
81 
93 
24 
38 
71 
26 
88 
23 
83 
09 
36 
90 
94 
63 
14 
56 
02 
57 
91 


Approximations to k, (s) 


Equation | 


(42) 


0.000 00 


0.108 
0.206 
0.297 
0.380 
0.456 
0.526 


0.590 3 


0.649 
0.703 
0.754 


0.800 3 


0.843 


0.882 


88 | 


67 


48 | 


0.919 
0.954 
0.985 
1.015 


1.043 3 


1.069 
1.093 
1.195 
1.273 
1.335 
1.387 


08 


93 | 


43 


| 
| 
| | 
|| 
| 
10 = | | 

09 | 
| 24 | = | 
-0.6 32 35 | 
—0.5 63 | 40 | 
~0.4 90 | 23 | | 
~0.3 | 78 — | 
-0.2 | 35 44 | 
-0.1 36 84 | 
0.0 | 3] | 02 | | 
2. 7 | 25 | 
0.4 4 | 69 80 | | 
0.5 61 04 | | 
0.6 31 | ZZ 04 | 89 | | 
0.7 55 | 23 | | 60 | | 
0.9 | 62 _ 29 | | 

1.0 | 62 | | 61 | 

1.5 | 2 | 39 | 

| 2.0 | 3 | BH 15 
25 | | 4 | | = 35 | 
3.0 | 16 | | 
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TABLE II (continued) 


4.0 


6.0 


7.5 


| 0.546 


0.600 
0.624 


0.725 


| 0.823 
| 0.869 
| 0.898 


U? (s) 
(Ref. 6) 


| 


0.436 
0.477 
0.514 


0.575 


0.645 
0.664 
0.681 
0.697 
0.712 


0.738 
0.760 
0.779 


0.916 
0.939 
0.953 
0.978 
0.995 898 | 
0.997 972 | 
1.000 000 


k, (s) 


| (Ref. 


1.436 


6) 


4 


1.477 7 


1.514 


1.546 
1.575 


1.600 
1.624 
1.645 
1.664 
1.681 
1.697 
1.712 
1.725 
1.738 
1.760 


1.779 
| 1.823 


1.869 


1.898 


1.916 
1.939 
1.953 


| 1.978 
| 1.995 


1.997 


2.000 000 


= 


972 


(s) 


P (s)=1 


for 


s>1 


Approximations to k, (s) 


Equation 


(38) 


1.425 84 
1.464 44 
1.498 16 

| 1.527 88 
1.554 24 
1.577 84 
1.599 04 
1.618 20 
1.635 64 
1.651 52 
1.666 12 
1.679 48 
1.691 88 
1.703 32 
1.723 80 
1.741 68 
1.783 64 
1.829 72 
1.859 60 
1.880 56 | 
1.908 04 
1.925 24 

| 1.961 36 
1.992 04 | 
1.996 01 

2.000 00 


Equation 


(41) 


1.437 
1.480 
1.518 
1.551 
1.579 
1.605 
1.627 
1.647 
1.665 
1.681 
1.695 
1.708 
1.720 
1.731 
1.751 
1.768 


1.810 3 


1.859 
1.895 
1.921 
1.956 
1.975 
1.998 


2.000 00 
2.000 00 | 
2.000 00 | 


00 
58 


65 


Equation 


(42) 


1.431 
1.471 


| 1.506 
| 1.539 


1.568 


1.596 


1.622 
1.646 
1.668 
1.689 
1.709 
1.727 


1.744 5 


1.760 
1.789 
1.815 
1.875 
1.934 
1.965 
1.982 
1.995 
1.998 
2.000 
2.000 
2.000 
2.000 


68 
00 
00 
00 
00 


| | | 

3.5 | 4 | Zs 78 
4.5 | 0 34 72 
55 | 0 96 | 94 
| = 8 23 s7 
65 58 26 
| 7.0 | 45 21 
| 21 59 
| 8.0 19 51 
| 8.5 | ri 64 10 
| 90 80 42 
| 10 oes 90 67 
| 12 89 48 
| 15 7 07 
20 | 88 78 
25 44 95 
| 30 | 81 23 
| 40 | a 23 16 
49 

100 | 00 

500 | 898 

00 
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8. Conclusion 


The aerodynamic theory, summarised in this paper is exact within the basic 
assumptions of two-dimensional incompressible aerofoil theory. Application of its 
results to the problem of gust loads on resilient aircraft has been given in a number 
of papers dealing with rigid and elastic aeroplanes. Because of the three-dimensiona] 
features of actual aircraft any such application can produce only approximate 
solutions. On the other hand it is known that the dynamic representation of 
structures of the complexity of aircraft is possible only in a very approximate 
manner. In view of this it appears that the results deduced here will be sufficiently 
accurate when dealing with the problem of gust loads on conventional low speed 
aircraft. This state of affairs will only be changed when some of the basic 
assumptions are no longer satisfied. For example, an investigation of swept-back 
wings cannot use these results since, quite apart from difficulties of structural 
representation, the flow around such wings cannot be approximated closely by 
two-dimensional conditions. 
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